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Theoretical and Experimental Studies Relevant to 
Interpretation of Auroral Emissions 


Introduction 


^ f increasin 8 awareness of the fragile balance required to 
maintain the Earths environment in a healthy condition. The Earth's upper atmosphere and its 

;“ Wlth f fK Un , ar f recognized as playing a crucial role in maintaining this balance 
Under tandmg of this link between the Earth and Sun contributes significantly to om 
understanding of Earth s atmosphere, climate, weather, magnetic storm activity, etc. Despite a 
rapidly increasing data base of knowledge on the solar-terrestrial environment much remafns to 
be learned Optical emissions from upper atmospheric constituents represent a fundamental tool 
for understanding both the composition and the physical processes that contribute to the ability 
of the atmosphere to sustain a high quality of life on the Earth. Space-based observation 
platforms are frequently able to make measurements of upper atmospheric emissions that are 
either impossjble to make from the ground or superior to their ground-based counterparts 
Specifically auroral imaging from a space-based platform makes it possible to determL the 
otol auroral energy influx and the characteristic energy of the incident auroral particles In 
addition it is possible through modeling, to map and relate these parameters from the 
thermosphere and ionosphere to the various regions of the Earth's magnetosphere 

This report describes the accomplishments of a program designed to develop the tools 
necessary to interpret auroral emissions measured from a space-based platform. The research 

ZerlnH T ^ “Ti ^ ** “ a laborator y Stud y d ^ncd to improve our 
derstanding of the space vehicle external environment and how it will affect the space-based 

measurement of auroral emissions. Facilities have been setup and measurements taken to 

simulate the gas phase environment around a space vehicle; the radiation environment 

encountered by an orbiting vehicle that passes through the Earth’s radiation belts; and the 

thermal environment of a vehicle in Earth orbit. The second major area of study was a modeling 

program to develop the capability of using auroral images at various wavelengths to infer the 

to al energy influx and characteristic energy of the incident auroral particles. An ab initio auroral 

calculation has been added to the extant ionospheric/ thermospheric global modeling capabilities 

wi hin our group. Once the addition of the code was complete, the combined model was used to 

comparethe relative intensities and behavior of various emission sources (dayglow aurora etc ) 

The work completed under this grant has provided an essential foundation for 

ImterTuvn th ( ™°\ a ' images that wil1 be available a ^r deployment of the Ultraviolet 
Imager (UVI) on the Global Geospace Science (GGS) POLAR satellite. Thus, a significant 

contribution has been made to the field of space-based auroral imaging. 

Cross Section Facility 

One of the ftmdamental factors that affects a space-based observing platform is the gas 
phase environment surrounding the platform. Atomic oxygen, in particular, plays a pivotal role 
in the induced environment for vehicles in low Earth orbit(LEO). Phenomena such as "shuttle 


g I° W ’" ,f 0sion of s P ace craft surfaces, and the “contaminant cloud” surrounding a spacecraft 
affect all instruments making observations from space. One of the tasks of this grant has been the 
development of a Cross Section Facility to address these issues. The Facility consists of a 
crossed molecular beam apparatus with a fast atomic oxygen beam and a thermal energy beam of 
another atmospheric gas (N 2 , 0 2 , etc.). A rotatable mass spectrometer serves as a sensitive 
particle detector The Facility is described in greater detail in the paper entitled “Laboratory 
Facility for Simulation of Vehicle-Environment Interactions” that was presented by the Principal 
nvestigator, Dr. Charles E. Keffer at the Vehicle-Environment Interactions Conference held at 
the Applied Physics Laboratory, Laurel, Maryland on March 1 1-13, 1991. The paper is attached 
as Appendix A of this report. 

... Init jfJ. Work with the Cross Section Facility focused on a thorough characterization of the 
acihty.. , Thls was essentiaI to ensure that measurements of gas-gas interactions such as 
differential scattering cross sections could be made without any serious systematic errors. 
Specifically, we have completed characterization of thermal energy beams of 0 2 and N 2 to 
determine their spatial and temporal characteristics. This information was used to assist in the 
measurement of the thermal energy differential scattering cross section for collisions of 0 2 on 
N 2 . Results of these measurements were compared with classical scattering theory and were 
found to be in good agreement (Figure 1). This serves as an indication that there are no serious 
systematic errors in differential scattering cross section measurements made with our Facility. 


°2 + N 2 



Center of Mass Angle 

Figure 1 : Differential scattering cross section for 0 2 on N 2 . 

Another major area of Cross Section Facility work was the development and 
characterization of an energetic, high flux source of atomic oxygen. The velocity and 
composition of the source were measured a large number of times. Beam velocity measurements 
are illustrated in Figure 2. Two peaks are evident in the intensity Vs time plot of the mass 
spectrometer signal. The first peak is due to photons from the plasma discharge formed from 
dissociation and ionization of molecular oxygen by a C0 2 laser pulse. These photons arrive at 



the detector essentially coincident with the formation of the fast O atoms. The second peak in 
the figure is due to the fast O atoms. The velocity of the atoms is the known distance from I the 
pulsed valve to the mass spectrometer divided by the time between the two peaks since this time 

" " me 1 fligh ' f ° r faSt ° a,0ms - Th * velocity for , he “to™ i 6 8 

km/sec. Mass spectrometer measurements of the fast oxygen beam demonstrated that it was 
composed predominantly of oxygen atoms. Some oxygen molecules, a small percentage of 
atomic oxygen ions, and some impurities are also present in the fast oxygen beL. Table l 
summarizes the characteristics of the energetic oxygen atom source. 
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Figure 2: Fast atomic oxygen pulse. 


Velocity Range 

3 to 1 1 km/sec 

Oxygen Atoms 

> 80 % 

Oxygen Molecules 

< 20 % 

Total Ions 

r< i % 

Impurities (H, N 2 , etc.) 

< i % ■ 


Table 1 : Energetic oxygen atom source operating characteristics. 

This work completes the development of a Cross Section Facility that is capable of 
performing a wide variety of vehicle-environment interaction measurements including both gas 
phase and gas-surface phenomena. This versatile tool should find extensive use in making 

aboratoiy measurements that can help in the interpretation of observations made from space- 
based platforms. F 


Induced Environment Workshop 

A workshop sponsored by this grant provided a forum for presentation and discussion of 


the latest research in the area of space vehicle induced external environments. A group of 
nationally recognized experts on space vehicle contamination met January 30-31 199? 
Huntsville, Alabama. Twenty invited participants each presented a paper in their particular area 
of expertise. The papers and discussion addressed three questions: 

(1) What is our current state of knowledge of the likely induced external environment for a 

large space-based platform in low Earth orbit? moment ror a 

(2) What progress has been made during the last two years in the vehicle contamination 

J^sTonsf ' “ * he predlctive capability of induced vehicle contamination for future 

smdies?™ 1 " iSSUeS remai “ unres0lvcd and " e the "><>* important to investigate in future 

Appendix B contains the summary report from the workshop, an agenda of the meeting 

of the attendees ' The Principal Investigator, Dr. Charles E. Keffer presented a paper 
entitled Cross Section Work at UAH/MSFC” at the workshop. P P P 

Radiation Testing 

As the Ultraviolet Imager (UVI) neared completion, the laboratory portion of this study 
began to concentrate on the effect of the natural and induced space environment on the UVI The 
UVI is a part of the International Solar Terrestrial Physics (ISTP) mission to study the Earth’s 
magnetosphere. It will be on a polar orbiting satellite exposed to the high radiation environment 
of space during its nominal three year mission lifetime. Of particular concern for the UVI is the 
all-dielectric multilayer interference reflection and transmission filters which are a part of the 
UVI optical system. The UVI filters will receive a radiation dose on orbit of < 250 krad 
Previous laboratory studies of high-energy radiation damage to optical materials at ultraviolet 

ToJa M gth i S 316 P T an y transmittance measurements of uncoated substrates (Heath and Sacher 
966 ; Nicoletta and Eubanks 1972; Reft et al., 1980; Pellicori * al, 1979; Hass and Hunter,’ 
1979 Gnllot and Rosenberg, 1989). Prior to the work under this grant, the extent of damage if 
any, that the UVI filters will likely experience from this radiation exposure was unknown We 

ave, therefore, performed a senes of tests to simulate the anticipated radiation dose and to 
measure its impact on the UVI filters. 

We selected 14 single layer thin films, 4 uncoated substrates, 2 multilayer transmission 
filters and 12 mult, layer reflection filters for radiatin damage study. A 60 Co source at Goddard 
S^ce fight Center irradiated seven of the single layer thin films, 2 of the uncoated substrates 
and all of the multilayer filters with 250 krad of gamma radiation. The remaining 9 samples’ 
served as controls with no irradiation. Table 2 shows a matrix of the 22 samples used with the 
substrate and thin film type shown together with the UVI filter to which it is similar The 

complete set of samples is analogous to the substrates, thin films, and filters that make up the 
UVI optical system. F 


Sample 

Substrate 

Coating 

Radiation/Control 

UVI Filter ] 

AI 2 O 3 1/2 

Fused Silica 

AI 2 O 3 

Radiation 

TSOLAR 

AI 2 O 3 2/2 

Fused Silica 

A1 2 0 3 

Control 

none 




Sample 

Substrate 

Coating 

Radiation/Control 

UVI Filter 

BaF 2 1/2 

MgF 2 

BaF 2 

Radiation 

T1304,T1356 TLBHS 

BaF 2 2/2 
CaF 2 1/2 
CaF 2 2/2 
F14-A 
Fused Silica 
uncoated 1/2 

MgF 2 

MgF 2 

BaF 2 

CaF 2 

Control 

Radiation 

none 

none 

MgF 2 

MgF 2 

Fused Silica 
Fused Silica 

CaF 2 

BaF 2 ,MgF 2 

Control 

Radiation 

none 

none 

none 

Radiation 

Tlbhl,tsolar 

r used Silica 
uncoated 2/2 

none 

Control 

none 

H 1 O 2 1/2 

MgF 2 

Hf0 2 

Radiation 

none 

HIO 2 2/2 

MgF 2 

Hf0 2 

Control 

none 

LaF 3 1/2 

MgF 2 

LaF 3 

Radiation 

~R1304,R1356,RLBHS,RLBHL, 
T1304,T1 356, TLBHS 

LaF 3 2/2 

MgF 2 

LaF 3 

Control 

none 

Mr45-C 

MgF 2 

BaF 2 ,MgF 2 

MgF 2 

Radiation 

none 

MgF 2 1/2 

MgF 2 

Radiation 

all 

Mgr 2 2/2 

MgF 2 

MgF 2 

Control 

none 

MgF 2 

uncoated 1/2 

MgF 2 

none 

Radiation 

none 

MgF 2 

uncoated 2/2 

P PI AD 

MgF 2 

none 

Control 

none 

-tx-T 1^3 .LV n 

RF135L/D 

Fused Silica 
Fused Silica 

LaF 3 /MgF 2 

Radiation 

none 

LaF 3 /MgF 2 

Radiation 

none 

Si0 2 1/2 ^ 

MgF 2 

Si0 2 

Radiation 

TLBHL 

01 O 2 2/2 

MgF 2 

Si0 2 

Control 

none 


Table 2. Summary of coatings and measurements. 

A senes of reflectance and transmittance measurements has been performed on the 22 
samples shown m Table 2. The preradiation measurements were performed in a hydrocLn- 
free cryopumped vacuum chamber at a pressure below 10' 5 Torr. A deuterium lamp with a MgF, 
window together with a 0.2-m vacuum monochromator provided 1 nm FWHM spectral 
resolution over the 120 to 1 80 nm wavelength range. Folding and collimating optics and a 6 mm 
diameter aperture limited the light incident on the eight position filter wheel holding the 12 5 mm 

appr ° ximately 1/4th * e of the thin film. A photomultiplier 
tube with a MgFj window and a semi-transparent Csl photocathode served as the detector for all 
measurements. The postradiation measurements were accomplished in another hydrocarbon-free 
ryopumped vacuum chamber also at a pressure below 10' 5 Torr. A similar deuterium lamp 0 2- 
m vacuum monochromator, and aperture were used for the postradiation measurements 
However the second chamber did not have any folding and collimating optics. The same Csl 
photomultiplier tube and eight position filter wheel was used for both the pre- and postradiation 



measurements. Reflectance measurements were performed at a 45° angle of incidence while 
transmittance was measured at normal incidence. Figures 3 and 4 illustrate r> f i 
configurations for the two vacuum chambers used. ^ 1Ca 

For the reflection measurements, two identical scans from 120 to 180 nm were performed 
on each substrate prior to their radiation exposure. Filter RF135L/G1I was measured aTa control 

me e RFi 1 35L«Ti eme T 7 check ' he . reliability and «P«>d»cibility of the data from scan 
'° , n - RF135L/GII control ensured that no systematic errors were introduced when the 

1 20 ml 7, ““ V , e "‘ Cd “ aew Set of subs,rates was installed. Two identical scans from 

FHte MnF n '**'* on each ™ b strate during the transmission measurements. 

Filter MgF 2 17 was measured on every transmission measurement as a check on the scan to scan 

reliability and reproducibility. The RF135L/GII and MgF 2 17 filters also provided traceabili* 
between the two vacuum chambers used in the pre- and postradiation measurements. ' 



Figure 4. Optical configuration for postradiation measurements. 

Repeated measurements of the RF135L/GII and MgF 2 17 filters have led to an estimated 
run to run uncertainty of approximately ± 3%. Of the 13 iiradiated samples, only 1 had any 

u7cMed M 7® T di f 1 rences in the 120 to 180 nm range after radiation exposure. The 
uncoated MgF 2 sample had an apparent 5% absoiption short of 1 60 nm before radiation exposure 
that was not present after radiation exposure. 




The primary wavelength range of interest to the UVI is 120 to 180 nm. However most 
ultraviolet optical materials that have been studied in the bulk form have their strongest 
ra lation-induced absorption bands at wavelengths longer than 180 nm. So another set of 
measurements was made from 160 to 300 nm where the irradiated samples w^e compared to 
their equivalent control samples because no preradiation measurements were made^n this 
wavelength range^A 3-5 /o radiation-induced absorption band centered at 260 nm was discovered 
m two of the MgF 2 substrates. This agrees with other published results on radiation damage to 
MgF 2 crystals. Since 260 nm is outside the wavelength range of interest to the UVI, however, 
this radiation damage is not a concern for the UVI. 

. In S ^ m f ry ’ lt appears that none of the samples tested have experienced any major 

the UVI 1 'slncffh nv e r S n, r ectance , 0r ,ransminance in wavelength range of interest to 
e UVI Since the UVI filters are made exclusively from materials that we have tested it is 

anticipated that they will suffer less than approximately 5% loss in reflectance or transmittance 
during their nominal 3 year mission lifetime. 

, Two pap< f rs were written on the radiation testing results for the UVI filter materials 

pr " S , ented 0n Jul y 16 ’ 1993 > in a paper entitled "Radiation Damage 

seSn of th^SP fT? Substrates " at the and Ultraviolet Polarimetry 

Dippn rJlT SPIE Inte 7 atlonal S y mposiu m on Optics, Imaging, and Instrumentation in San 
D^go, California. Complete results have been presented in a paper entitled "Radiation Damage 
Effects in Far Ultraviolet Filters, Thin Films, and Substrates" that is scheduled for publication in 

Applied Optics in 1994. Copies of these two papers are attached to this report as Appendix C 
and Appendix D, respectively. FF 

Thermal Testing 


to the t o, U aV10, t T er 15 paSS ‘ Vely C00led by rad ‘ at ' ve heat transfer from the instrument 
to the 3 K thermal background of space. The performance of the UVI in this thermal 

environment i> a key factor in the success of the mission. The radiative coupling between the 
UVI and the 3 K thermal background of space must meet the critical cooling^equirements for 
the instrument. These include a CCD detector temperature below - 55 °C for noise reduction an 
optical bench temperature near room temperature to prevent misalignments caused by thermal 
expansion or contraction, and heat removal from "hot" electronic parts to prevent failures We 
have, therefore, designed and earned out several tests of the UVI thermal performance. 

A test chamber setup that was designed to simulate the thermal environment of the UVI 
on orbit was built into the VUV calibration vacuum system. The UVI was mounted on a 
baseplate cooled or heated over the nominal temperature range of - 20 to 30 °C by a recirculating 
water and ethanol mixture. A second cold plate was suspended approximately 2" above the UVI 
ra rntor. This cold plate was painted black to increase its emissivity and thus improve the 
radiative heat transfer from the UVI to the cold plate. Liquid nitrogen was flowed though the 
upper cold plate so that its temperature approached 77 °K. Several thermistors and 
thermocouples monitored the temperature at strategic locations. 

Several key questions were addressed in this thermal test. Will the CCD reach its 
operational target temperature of < - 55 °C? Can the optical bench be maintained near room 



temperature with 10 watts of available heater power? Will the electronics remain sufficiently 
warm dunng survival mode conditions? Do the "hot" electronic parts have adequate cooling? In 
cases, the UVI thermal performance was found to be satisfactory. All the required thermal 
vacuum cycles and thermal balance tests were also completed. 

UVI Optical System Testing 

at f The 3 JVI IS desi § ned t0 measure atomic oxygen auroral emissions at 1 30.4 and 1 35 6 nm 
N 2 Lyman-Birge-Hopfield (LBH) auroral emissions from 140 to 160 nm and from 160 to 180 
nm, and the scattered sunlight from 180 to 210 nm. High sensitivity at the selected wavelength 
coupled with a high rejection ratio for other wavelengths is required. The UVI uses a filter wheel 
with five filter sets for the five different wavelength ranges. Each filter set consists of 3 
reflection filters and 1 transmission filter. 

The success of the UVI optical system depends on the correct design and performance of 
these filters. Accordingly, measurements have been completed of the reflectance or 
transmittance of all of the flight filters to verify that they achieve the necessary performance 
requirements. The measurements were completed in the same manner as the radiation test 

measurements and using a setup similar to that illustrated in Figure 4. The net transmittance for 
each filter set is shown in Figure 5. 



Figure 5: UVI flight filter sets. 

In addition to the five filter sets, the UVI optical system has three telescope mirrors a 
MgF 2 window in the aperture door, and a folding mirror for diverting light from the primary 
detector system to a secondary detector system if desired. Reflectance measurements have been 
completed on test coupons coated at the same time as the folding mirror and each telescope 
mirror. Reflectance values of 70% are typical in the 120 to 210 nm range. Transmittance 
measurements have also been completed on several conductive films designed for use on the 
MgF 2 window in the aperture door. 


Auroral Modeling 


An unambiguous interpretation of auroral emissions is dependent on a successful 
modelmg program. In particular, the modeling program must explicitly include the arnora 
electrons and the assoc, ated auroral emissions while maintaining the capability of Lcorpomdng 
the auroral results within a comprehensive global ionospheric model. The long term goal of om 
modelmg program ,s to develop the capability of using auroral images at various wavelet Asm 
infer the total energy influx and characteristic energy of the incident auroral ^ lef ,n 

“;ri: e corc,':ct able t0 dCdUCe — -position and 

• In 0 J“" research plan we proposed a program to incorporate an auroral code into our 
wistmg global model and conduct studies with the resultant computer model. As an associated 
part of the second task, we also provided support for two meetings of the members of the ISTP 
Ultraviolet Imager (UVI) science team. As described below, we have successfully met these 

Integration of Models 

We began this program with a global ionospheric model (FLIP) in a fullv 
interhemtsphenc format, including photoelectron transport. The model accounted for the 
chemistry of all major emitting species, solved for key mechanisms, and provided the hill odd 

1990;Toit ' 1983; Torr and Terr, 1982; Younger 
, 1980). The global model, however, did not include the capabilities to model auroral 
e ectrons and their resultant emissions. The treatment of photoelectron transport within the 
global model does, however, readily lend itself to inclusion of auroral electrons. In preparation 
of our long torn modelmg needs, we had previously developed (with P. G. Richards) a two- 
stream auroral electron energy code that determines the energy degradation of the primary 
spectrum as a ftmc “ on °f energy and altitude; it is fully described in Richards and Torr, r 1 990 1 
The integration of our two-stream auroral code wi A the global FLIP model is central to 
Ae modeling aspect of A,s program. The auroral electron code accurately models energy 
degradation of an incident electron flux and resultant auroral emissions. It has no knowledge 
however, of the chemistry and dynamics of the global ionosphere-thetmosphere. The FLIP is m 
excellent representation of Aese processes but does not include the auroral processes The 
combined model allows simultaneous modeling of auroral emissions superimposed on 
un er ytng airglow. The updated FLIP model now accepts auroral production data from Ae 

capabditiies de “ Standard mPU ‘ ^ “ d ‘ S available for studies ‘hat exploit its combined 

High Latitude Extensions 


The FLIP model 1S a closed field line model and assumes, incorrectly, that the field lines 
remain closed throughout the polar regions. Attempts to model these regions with closed field 



ines lead to infinitely large field lines and attendant computational difficulties We had 
previously thought the FLIP model would have to be modified to include open field lines Such 
a modification would have been an extensive project. Our work under this grant has shown 
however, that a good model of the airglow can be obtained by turning off the field-line diffusive 

modification hC ^ ? ^ ^ ^ ^ advantage of using the existin 8 code with only minimal 

In addition to being able to run at polar latitudes, the FLIP model must also be modified 
to take account of several high latitude mechanisms. One of the more significant mechanisms is 
the drift of ionospheric plasma due to local electric fields. Plasma convection throughout the 
polar region can significantly alter the plasma densities relative to a static model which has the 
feet of the flux tubes anchored to the Earth’s surface. Plasma convection is controlled bv 
ionospheric electric fields which are intimately connected with the magnetospheric-ionospheric 
current systems and are generally not readily modeled because of their complexity. There are a 
number of empirical algorithms, however, which accurately model the overall convection 
patterns. We adapted the convection model of Sojka er al. [1986] as a first step to allow the 
FLIP model to include plasma drift. We are currently able to specify plasma trajectories as a 
function of time and interplanetary magnetic field strength and direction. This capability is 
significant in its own right and can be used in numerous studies. It will also form the basis for 
future extensions to the FLIP code. 


Modeling Studies 


... . The corn ^ med codes h ave enabled studies of superimposed airglow and aurora. In 
addition, the models have been used to calculate conductances due to both EUV and auroral 
sources. Details of these studies are given below. 

FUV Auroral Emissions 


Rees and Luckey [1974] showed that measurements of the 6300 A emission of atomic 
oxygen and the 4278 A emission of ionized molecular nitrogen could be used to infer the 
characteristic energy of the incident auroral particle flux. In this technique, the energy 
dependence is provided by the collision^ deactivation of the 6300 A emission at lower altitudes 
Strickland et al. [1983] demonstrated the same capabilities using FUV auroral emissions In 
particular, they proposed the use of the O I (1356 A) and N 2 Lyman-Birge-Hopfield (LBH) 
emissions as remote diagnostics of the auroral process. Instead of quenching mechanisms, the 
FUV technique relied on absorption of 0 2 at lower altitudes to supply an energy dependence. In 
a previous study [Germany et al., 1990], carried out under contract NAS8-37586, we examined 
the stability of FUV emission ratios formed from atomic oxygen and LBH emissions as a 
function of likely changes in the neutral atmosphere. These are the FUV emissions we studied 
under this grant and are shown in Figure 6. Some emissions (O I 1356 and LBH 1464) vary 
significantly with average energy, while others (LBH 1838) are relatively independent of energy 
an are indicators of the total flux. They found that emission ratios formed from the O I (1356 
A) emission varied significantly with seasonal and solar cyclic variations. On the other hand, 
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Figure 6. FUV auroral emissions used as 
diagnostics of the auroral processes. 
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Figure 7. Ratio of FUV auroral emissions as a function of average 
energy. The 1356 A emission ratio has a greater sensitivity, but also 
has a greater variability than the LBH ratio. Dashed lines ate for 
summer conditions (day 173) and solid lines for winter (day 356). 
Fot each day, solar activity is varied by changing the F107 and Ap 
indices. The largest intensity ratios are obtained at solar minimum 
activity; the smallest ratios occur at solar ma-rim™™ 




emission ratios employing long- and short-wavelength LBH emissions were much less sensitive 

LBHmSl'i?* neUtral atm0Sphere (Figure 7 >' Furt Itermore, blending of the O I (1356 A) and 
LBH (1354 A) emissions represents a serious observational difficulty. Consequently the LBH 

emission ratios are potentially more valuable as remote diagnostics of the incident auroral flux 
than similar ratios using the O I (1356 A) emission. 

Much of our work under this grant on the use of auroral LBH and 01 emissions as remote 
diagnostics of auroral processes is summarized by Germany el al. [1994], A reprint of this paper 

' Veexan,ined the dependence of these emissions and their ratios on 
e shape of the energy distribution of the auroral electrons. We were particularly concerned that 
the column-integrated emissions might be overly sensitive to the energy distribution of the 
mci ent energy flux. Since different physical processes lead to different energy distributions 
there is not a single parameterization which applies to all situations, furthermore since the 
energy distribution is not typically known for remote sensing applications, a dependence on the 
d stnbution could severely limit the information which could be retrieved from the emissions 
Strickland et al. [1983] examined the dependence of FUV auroral emissions on the choice of 
either a Gaussian or a Maxwellian incident energy spectrum. We extended that study by 
modeling auroral emissions resulting from incident energy distributions of different widths and 
shapes to determine whether changes in energy spectra might obscure, or render ambiguous, our 
determination of the incident characteristic energy from FUV emission ratios. In addition to 
varying the functional dependence of the energy distribution, we looked at the effects of varying 
its width and its low energy behavior. Results are shown in Figures 8-10. 

Fortunately, we found that modeled column-integrated emissions showed less than 30% 
variance with changes in the shape and width of the incident energy spectrum, provided the 
average energy and total energy flux of the energy distribution are held constant. Our choice of 
istributions was not necessarily realistic; the distributions were chosen to illustrate extremes of 
shape and width. It is significant, in light of these choices, that a change of energy distribution 
affects the column emissions by less than 30%. This is relatively small, given the other 
uncertainties associated with auroral imaging. Regardless of the energy distributions modeled, 
we found that the differences in resulting emissions are governed primarily by the average 
energy of the distributions and the relative widths of the distributions. This is principally duelo 
the predominate influence of 0 2 absorption in determining the energy dependence of the FUV 
aurora emission. Changes in ratios of FUV emissions as a function of characteristic energy are 
much larger than those due to changes in energy distribution. Consequently, the uncertainties 
associated with the energy distribution are not large enough to significantly affect the 
characteristic energy determination. 

The combined auroral and FLIP models allow simultaneous modeling of the EUV and 
auroral souce functions We modeled typical aurora with concurrent dayglow emissions as a 
function of solar zenith angle, solar activity levels, total electron energy flux, and incident 
charactenstm energy. Figure 1 1 shows modeled auroral and dayglow emissions for O I (1356 A) 
and LBH (1464 A). The ratio of auroral to dayglow emissions was found to be a strong function 
of solar zenith angle because of the dayglow dependence on solar zenith angle. We found that 
the airglow can be quite intense relative to auroral emissions. Indeed, under certain 
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Figure 8. The top panel gives die ratio of column 
emissions from two energy distributions: a Gaussian 
and an Maxwellian (bottom panel). Both energy 
distributions have the same average energy and total 
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Figure 9. Same as previous figure, except for 
Gaussian distributions of varying width. Both 
energy distributions have the same average energy 
and total flux. 






Figure 10. Same as previous two figures except for 
two arbitrary energy distributions. Both energy 
distributions have the same average energy and total 
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Figure 11. Simultaneously modeled auroral 
emissions with underlying airglow. Solar activity is 

varied by changing the F107 and Ap indices. 





condmons (energy flux - 1 erg cm’ 2 s' 1 and solar zenith angle = 50 degrees) the dayglow will be 
the limiting factor in the detection of weak auroras. 

Auroral Conductances 

The combination of the FLIP and auroral codes also allows the calculation of ionospheric 
conductances due to auroral precipation. Knowledge of ionospheric conductivities is central to 
an understanding of the coupled nature of the ionosphere, thermosphere, and magnetosphere. 
This is reflected in models of ionospheric electrodynamics and convection requiring conductance 

from FlT , r derable eff0rt haS been direCted t0Ward “e i°nospheric^conductances 
from EUV solar illumination, nocturnal sources, and auroral precipitation. 

A key t0 remote determination of auroral conductances by optical methods is the ability 
to characterize the incident auroral flux. Since we have demonstrated the resultant colurn^ 
intensities are relatively insensitive to choice of energy distribution [Germany et al. 19941 the 
principal task is to estimate the energy flux and characteristic energy. 

In A ltial '' esults of our conductance modeling were presented at the 1 99 1 Fall meeting of the 
^ U . An abstra * of that Presentation is included in Appendix F. In addition, we have 
submitted a paper for publication discussing our use of LBH emissions to determine auroral 
conductances^ This paper is currently in review; a preprint is included in Appendix G The 
modeling study presented in the paper is a continuation of our previous examinations of the 
u i ity of LBH emission ratios as remote diagnostics of auroral processes. In this work, we 
examined the viability of using LBH emission ratios to infer auroral conductances and to 
quantify the strengths and weaknesses of this technique. To our knowledge, this is the first time 
e LBH emission ratio has been used to infer auroral conductances. We demonstrated that 
determinations of average energy and energy flux from a single measurement of two LBH 
emissions may be used to estimate column-integrated ionospheric conductances. Figure 12 
shows modeled ionospheric conductances as a function of the ratio of two LBH emissions We 
have examined changes in activity levels by varying the FI 07 and Ap indices in our models. 

Conductances are directly proportional to the corresponding ionospheric electron 
densities. We calculate conductances using the densities, temperatures, and magnetic field 
strengths provided by the FLIP model. The dependence of the determined conductivities on 
incident average energy, total energy flux, and changes in solar and magnetic activity levels was 
examined. We showed that, for energies above 5 keV, auroral conductances may be scaled bv 
he square root of the incident energy flux with errors less than about 20%. For lower energies 
however the scaling may deviate significantly from a power of 0.5. We provided appropriate 
scaling factors as a function of average energy. We also noted that the choice of either a 
Gaussian or Maxwellian distribution can significantly affect low-energy (<1 keV) conductance 
calculations (Figure 13). 

Finally, we quantified the conditions under which the mean energy of a Gaussian energy 
istnbution may differ significantly from the characteristic energy. It is common to state that the 
mean energy of a Gaussian energy distribution is its central, or characteristic, energy Since 
energy is constrained to be non-negative, however, the distribution may be truncated at zero 
energy. If so, the distribution will be assymetric and the mean energy will be larger than the 
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Figure 12. Ionospheric conductances as a 
function of the FUV LBH emission ratio. Solar 
activity is varied by changing the F107 and Ap 
indices. 



Figure 13. The choice of energy distribution can 
affect the modeled ionospehric conductivities. 

Note that the conductivity due to nighttime sources 
(starlight and reflected EUV) can exceed that from 
auroral sources. 





characteristic energy. The approximation of the mean energy by the characteristic energy fails 
or energy distributions that are very wide compared with the central energy. The quantitative 
meaning of the phrase ‘very wide' is given in the preprint contained in Appendix G. 

1356 A from 0 2 Dissociative Excitation 


In addition to the combination of the auroral and FLIP models, we’ve modified the 
auroral code to include an additional source of auroral 1356 A emission. Most models, including 
our previous codes, have only considered a single source of 0(5S) 1356 A emission namely 
electron impact excitation of O. There is, however, an additional potential source of 1356 A 
emission, the dissociative excitation of molecular oxygen by electron impact. Dissociative 
excitation of 0 2 by electrons is a significant source of O I 7774 A emission for hard auroral 
spectra [Strickland, et ai, 1989], We could reasonably assume that the same mechanism may 
also be important in the production of auroral 1356 A emission. Dissociative excitation cross 

sections for 1356 A have been reported [Wells et ai, 1971] but these have not generally been 
included m auroral codes. J 


Based on our initial modeling, we estimated factor of two changes in auroral 1356 A 
column intensities due to this additional source. Our initial work, however, failed to account for 
su suquent revisions of calibration standards. As a consequnce of these revisions, the 1356 A 
column enhancements due to the additional excitation source are about 15% at 20 keV incident 
energy. The cross section picture is somehat unclear and is detailed below. 

As noted above, the cross section for the production of 1356 A emission by electron 
impact of 0 2 has been revised since its initial estimate. Ajello [1971] measured the relative cross 
section lor 1356 A emission from electron excitation of 0 2 , estimating a lower limit of 1.65 x 10‘ 
cm for the peak cross section at 100 eV. Central to this determination was the assumption of 
thermal velocity distributions for the dissociation fragments and the use of the radiative lifetime 
of the 0(5S) state (565 ps) reported by Garstang [1961], Freund [1971] measured relative 1356 
A emission functions and, along with Borst and Zipf [1971], observed non-thermal velocity 
distnbunons for 0(5S) by time-of-flight (TOF) methods. Based on these studies, Wells e/ ai 
11971] subsequently used a radiative field-of-view analysis to combine the work of Ajello 
[1971] and Freund [1971] and estimated a peak value of 7 x 10' 18 cm 2 ± 50% at 100 eV They 
then renormalized the Ajello [1971] and Freund [1971] relative cross sections to this value. 

The radiative lifetime of the 0(5S) state was later measured and revised downward by 
about a factor of 3 [Johnson, 1972; Wells and Zipf, 1972; Wells and Zipf, 1974], The 1356 A 
peak cross section should therefore be revised downward by the same amount. Using the latest 

value of Wells and Zipf [1974] we estimate the peak cross section at 100 eV to be 2 1 x 10‘ 18 
cm . 


We note that none of the lifetime determinations [Johnson, 1972; Wells and Zipf 1972- 
Wells and Zipf, 1974], the TOF studies [Freund, 1971; Borst and Zipf, 1971], or the analysis of 
Wells et ai [1971] were dependent on the Mumma and Zipf [1971] Lyman-alpha calibration 
benchmark. Ajello [1971] did employ a similar calibration, but it was independent of that of 
Mumma and Zipf [1971], Hence subsequent revision of this benchmark [Shemansky et ai. 



1985; Van Zyl et al., 1985; Ajello et al., 1988; Ajello et al., 1991] should not have affected the 
1356 A dissociative emission cross section. 

There appear to have been no recent measurements of the cross section. Though we have 
found no references to additional work beyond that discussed above, there are some conflicting 
references. Both Itakawa et al. [1989] and Ajello and Franklin [1985] reference the work of 
Wells et al. [1971] without mentioning the later lifetime revisions. On the other hand, Erdman 
and Zipf [1987] reference the lifetime revisions and then suggest an additional revision due to 
changes in the Lyman-alpha calibration benchmark. As noted above, we feel this is not 
necessary. 

Clearly the 1356 A dissociative cross section could benefit from additional study. 
Regardless of the value used, any dissociative emission cross section that is an appreciable 
fraction of the atomic cross section may contribute substantial 1356 A emission due to the broad 
nature of its excitation function. This is the basis for our modeling of the effects of this 
additional source, using the cross sections of Wells et al. [1971] revised to reflect the lifetime 
determination of Wells and Zipf [1974], 

UVI Workshops 

The modeling task included support for an Ultraviolet Imager (UVI) science working 
group meeting in the second and third years of support to enable prelaunch planning and to gauge 
science team requirements. It readily became apparent, however, that there were several issues 
that needed immediate discussion by the group. Therefore, the first UVI workshop was held on 
15 & 16 August 1990 in Park City, Utah. The attendees included nationally prominent scientists 
in the fields of ionospheric and magnetospheric physics and chemistry and are listed in Table 3 
along with the meeting’s agenda. A major focus of the workshop was the mission science 
planning and the specific science requirements of each participant for the UVI instrument. The 
workshop was deliberately informal, with each participant giving presentations about topics of 
interest to them. 


A 

genda 

August 15, 1990 

August 16, 1990 

Introduction 

Discussion of Analysis Tools 

Instrument Observing Sequences 

Discussion of Signal Extraction 

Review of Science Objectives 

Revised Data Analysis Plan 


Summary 


Attendees 

Dr. Joe Ajello 

Jet Propulsion Laboratory 

Dr. Ken Clark 

University of Washington 

Dr. Bob Clauer 

University of Michigan 

Dr. Glynn Germany 

University of Alabama in Huntsville 

Dr. George Parks 

University of Washington 

Dr. Jim Spann 

^NASA/MSFC 




Dr. Doug Torr 

University of Alabama in Huntsville 

Dr. Marsha Torr 

NASA/MSFC 


Table 3. Attendees and agenda of the 1990 UVI workshop held in Park City, Utah. 

The second Ultraviolet Imager (UVI) workshop was held on August 1,2 1991 in Salt 
Lake City, Utah. The second workshop built on the progress of the first meeting and allowed 
each of the attendees to participate in the mission science planning and to address their specific 
science requirements for the UVI instrument. A major focus of the workshop was a review of 
available thermospheric/ionospheric modeling and image analysis techniques. The second 
workshop participants and agenda are given in Table 4. 


Agenda 

August 1, 1991 

August 2, 1991 

Goals for Auroral Imaging 

Program Status and Introduction 

VUV Calibration Source 

Automated Auroral Image Analysis 

Improved Neutral Wind Modeling 

Extraction of Parameters from Images 

Summary 

Update on Data Analysis 


Auroral Measurements and Models 


Attendees 

Dr. Joe Ajello 

Jet Propulsion Laboratory 

Dr. Keith Cole 

La Trobe University, Australia 

Dr. Ken Clark 

University of Washington 

Dr. Glynn Germany 

University of Alabama in Huntsville 

Dr. Phil Richards 

University of Alabama in Huntsville 

Dr. Ramin Samadani 

Stanford University 

Dr. Stan Solomon 

University of Colorado 

Dr. Doug Torr 

University of Alabama in Huntsville 

Dr. Marsha Torr 

NASA/MSFC 


Table 4. Attendees and agenda of the 1991 UVI workshop in Salt Lake City, Utah. 


The planning and coordination of both meetings was the responsibility of Dr. Germany. 
His responsibilities included arrangement of meeting facilities, communication with the science 
team members, and preparation of post workshop mailings. In addition. Dr. Germany 
participated in the workshop discussions and gave presentations at both meetings. 
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ABSTRACT 

A facility for simulation and study of interactions between a spacecraft in low 
Earth orbit (LEO) and its ambient environment is described. The facility is composed of 
a crossed beam apparatus with a rotatable mass spectrometer detector. It can be used for 
a wide variety of vehicle interaction studies including both gas phase and gas-surface 
interactions. Measurements of differential scattering cross sections, surface scattering 
phenomena, and spacecraft glow are representative examples. A key element in the 
facility is a laser-induced discharge energetic oxygen atom source for simulation of the 
ambient vehicle environment. Measurements of important characteristics of the oxygen 
atom source, including velocity and beam composition, are presented. Performance of 
differential scattering cross section measurements is evaluated using low angle scattering 
from thermal energy collisions between beams of oxygen and nitrogen molecules. 


INTRODUCTION 

Spacecraft in low Earth orbit (200-700 km) are exposed to an intense flux of 
atomic and molecular species from the atmosphere. At these altitudes the major 
constituent of the atmosphere is atomic oxygen which is created by solar UV 
photodissociation of molecular oxygen. Typical orbital velocities of 8 km/sec produce 
an atomic oxygen kinetic energy relative to the spacecraft of nominally 5 eV. The 
density of atomic oxygen varies with the amount of solar activity, but an average value 
for 250 km is about 10 9 atoms/cm 3 The flux of atomic oxygen impinging on spacecraft 
surfaces under these conditions is thus approximately 10 15 atoms/cm 2 /sec. This 
environment is known to cause a variety of phenomena including severe surface erosion 
of polymers and some metals as well as induced optical emissions in the ultraviolet, 
visible, and infrared portions of the electromagnetic spectrum [Bareiss et al., 1987], 


Since the discovery of these effects several years ago, a significant effort has been 
made to study them. A combination of spacecraft glow data (see, e.g., Torr et al., 1977; 
Banks et al., 1983; Mende et al., 1988] and several recent laboratory investigations 
[Arnold and Coleman, 1988; Caledonia et al., 1990; Holtzclaw et al., 1990; Orient et al.. 
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Dhenom^ a led H° ** develop ™ eilt of a preliminary data base for the spacecraft glow 
P a. However, very httle is yet known about such fundamental atomic and 

molecular parameters as scattering, excitation, or ionization cross sections for 5 eV 
ygen colliding with other species. Mechanisms for the spacecraft glow phenomena and 
the surface chemistry processes [Kofsky and Barren 1986] are just beginning to^ 
understood with much still remaining to be learned. ® 

An improved understanding of these atomic oxygen effects on orbiting spacecraft 
is important to the success of Space Station Freedom and to future Soace Shurflf* mice 
w* W*- which arc sens.nve ,o the indued cnwolcm^rh vchid ~ 
oxygen also plays a fundamental role in the chenusuy of the upper aimSere aXn 
high temperature combustion reactions. " 


FACILITY DESCRIPTION 

Figure 1 shows a schematic diagram of the facility. The basic design consists of 

which ’ ° g “ bCa T sources 311(1 311 in *plane quadrupole mass spectrometer detector 
which is rotatable around the collision region of the two gas beams [Lee et al 1969] 

The beam sources are pulsed to minimize pumping requirements and to optintize the 
energetic oxygen atom source in one of the beams. The configuration shown illustrates 
e setup for making differential scattering cross section measurements For surface 
scattenng measurements only one beam valve is used and the surface materid is placed 
at the center of rotation of the mass spectrometer detector. 

Vacuum System 

The cylindrical vacuum chamber is fabricated from 304L stainless steel with an 
inside diameter of 47" and height of 25.5" (see Figure 2). The walls Zmade thick To 
unit distortion and misalignment due to atmospheric pressure. The lid which is 

“ le . l ° aCCCSS t0 the „ interior of toe chamber, is 1.5" thick, the bottom is 1" 
211 e cylinder wall is 0.5" thick. The main chamber is pumped by a 3000 1/sec 
cryopump which removes essentially all of the gas from each pulse prior to 
of the next pulse at the usud 0.5 Hz repetition rate of the waives. ? “ 

, , „ H 1 ' ? T 5C£ f VaIVCS ^ in d * ffirendaI1 y pumped chambers each evacuated by a 1000 

1/sec turbomolecular pump. A 1 mm diameter beam skimmer separates each pulsed vdve 

chamber from the main chamber and coUimates the gas beam* Flexible bellows inlhe 

tiie S aa V 1 VC Chamb ™’ ! ogether with horizontd and vertied adjusting assemblies dlow 
the gas beams to be aligned so as to collide at the center of rouition of Ae TZ 

spectrometer detector. Precision beanngs and shafts in the adjusting assemblies dlow an 
ahgnmem accuracy of ± 005". Assure in the pu.sed valve cham4 “ 

Torr dunng operanon to limit loss in beam intensity due to scattering. 

nnmn “f™ ] S / pectrometer is ™unted inside of a housing which is differentially 
P ped by a 230 1/sec ion pump. Interchangeable apertures on the front of the housing 
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Figure 2. Sectional view of vacuum chamber. 
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and the entrance to the mass spectrometer ionizer limit the acceptance angle of the mass 
spectrometer to the region of overlap of the two gas beams. The mass spectrometer 
housing is suspended from a rotary vacuum seal which is differentially pumped by a 
small mechanical pump and a small ion pump. This arrangement allows the mass 
spectrometer to be rotated about the scattering center through approximately 200 degrees 
without significantly effecting the pressure in the main chamber or in the mass 
spectrometer housing. Base pressure in the mass spectrometer housing is 10' 9 Torr. No 

measurable effect on the pressure in the mass spectrometer occurs from operation of the 
pulsed valves. 


Atomic Oxygen Source 

Atomic oxygen is formed in one of the beam sources by pulsed laser-induced 
breakdown of pure molecular oxygen in the manner developed and described by 
Caledonia et al 1987. A sectional view of the atomic oxygen source is shown in Figure 
3. A pulsed valve in the center of the figure introduces a short pulse of molecular oxygen 
into a conical expansion nozzle. An approximately 8 J C0 2 laser beam passes through 
a 1 meter focal length BaF 2 lens shown at the lower left of the figure. The CO laser 
beam then passes through an AR coated ZnSe window which separates the vacuunAn the 
pulsed valve chamber from atmospheric pressure. The beam is focussed to about 1 cm 2 
on a gold coated nickel copper mirror. A smaller spot size results in an energy density 
on the mirror which can cause significant damage to the mirror coating. The 50 cm 
radius of curvature of the gold coated mirror focusses the C0 2 laser beam to about 1 mm 2 
near the orifice of the valve. A direct hit on the valve orifice is avoided since this 
damages the o-nng tip which seals the valve closed between pulses. The energy density 
in the focussed laser beam is sufficient to dissociate the molecular oxygen creating a high 
temperature plasma near the throat of the nozzle. The hot plasma accelerates down the 
expansion nozzle as a blast wave dissociating and accelerating most of the remaining 
molecular oxygen. Fast oxygen atoms exiting the nozzle are collimated by a 1 mm 
diameter skimmer as they enter the main vacuum chamber. Residual gas deflected by the 
skimmer is removed from the pulsed valve chamber by the turbomolecular pump. A 
small fraction of the fast gas pulse is composed of ions. These are readily removed from 
the pulse by deflection plates in the main vacuum chamber. These deflection plates are 
spaced 1 cm apart and biased with a few hundred volts poten tial difference. 


FACILITY PERFORMANCE 


Differential Scattering 

Classical scattering theory predicts that an intermolecular potential of the form 

V(r) oc -l/r* (i), 

where r is the intermolecular distance, should have a low-angle differential scattering 
cross section equal to 6 
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1 ( 0 ) oc 0 -™ 


( 2 ), 


where 0 is the ccnter-of-ma ss scattering angle [Bernstein, 1964], This functional form 
of the potential is typical of the long range attractive van der Waals forces between atoms 
and molecules. Experimental results presented as log-log plots of differential scattering 
cross section vs center-of-mass angle have confirmed this predicted behavior in the range 
of center-of-mass angles from roughly 1 to 10 degrees for a large number of collision 
partners [Bernstein, 1964]. We have used this simple relationship to validate the 
capability of our facility for making differential scattering cross section measurements 
without senous systematic errors. Collisions between beams of thermal energy oxygen 
and nitrogen molecules have been used as a test system. Several sets of data were taken 
or center-of-mass angles less than 10 degrees. Transformation of the data from the 
laboratory to the center of mass coordinate system was accomplished assuming the 
approximation of elastic scattering of monoenergetic beams colliding at 90° with no 
divergence and scattering detected in the plane of the two beams [Bernstein, 1964], 
Average values for the relative differential scattering cross section are shown in Table l' 


Table 1. Relative differential scattering cross section for 0 2 + N 2 collisions. 


Center of Mass Angle 

Relative Differential Scattering 


Cross Section 

(degrees) 

(arbitrary units) 

6.4 ± 0.2 

26.58 ± 0.42 

7.0 ± 0.2 

19.91 ± 0.08 

• 7.5 ± 0.2 

17.70 ± 0.42 

8.0 ± 0.2 

14.13 ±0.27 

8.5 ± 0.2 

13.21 ± 0.36 

9.1 ± 0.2 

11.75 ±0.14 

9.6 ± 0.2 

10.03 ±0.17 


Uncertainties in the angles are due primarily to the positional accuracy of the rotary 
vacuum seal. . Relative differential scattering cross section uncertainties represent 1 
stan d deviation for each set of data at each angle. A least squares fit to the log of 
the scattered intensity vs the log of the center-of-mass angle results in an experimental 
value for the exponent of -2.30 ± .13. This is within one standard deviation of the 
theoretical value and so indicates satisfactory agreement. 
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Atomic Oxygen Source 


Measurements have been made on the atomic oxygen source to determine the 

^ at0m , beam ' FlgUre 4 indicates type of measurements which have been 
made. The figure shows three peaks. The first one is due to photons from the oxygen 
plasma discharge which strike the Channeltron detector of the mass spectrometer while 
““ second peak is due to fas. oxygen atoms. The third peak on thHhoXr „f ti!e 

?“* - S dU , C Sma11 amount of Preamp ringing arising from the large, fast 
se time O atom signal. The photon peak serves as a time marker since the arrival of 
e photons at the detector is essentially coincident in time with the formation of the fast 
atoms. So, the time between the two peaks represents the time of flight for the oxygen 
atoms The velocity is thus calculated from the time of flight and the known distance 
from the valve to mass spectrometer. The average velocity is 6.8 km/sec for the oxygen 
atom pulse shown in Figure 4. The velocity of the oxygen atoms can be varied by 
changing the laser intensity or the amount of gas per pulse. A velocity range of 3 to 1 1 
knvsec has been measured for this source. 

The fast atom beam composition has also been measured to determine the ratio 
of oxygen atoms to oxygen molecules and to identify the presence of any impurities in 
bea ?'oA M ff S trometer measurements of the beam have been made indicating a 
range of 80 % to 95 % atoms in the beam. Two competing processes in the mass 
spectrometer ionizer contribute to an uncertainty in the actual fast atom content of the 
beam as it is deduced from the measurements. Dissociation of molecular oxygen by the 
ionizer results in an overestimation of the atom content in the beam. Ten percent 
dissociation of molecular oxygen is typical in our ionizer. For a beam which is 

0Xygen ’ *“* rcpresems a smaU uncertainty. For example, a beam 
of 80% O and 20% 0 2 would appear as 82% O and 18% 0 2 in the mass spectrometer due 
to e dissociative ionization. On the other hand, recombination of atoms in the mass 
spectrometer during the duration of the fast pulse results in a loss of atoms in the beam. 
The degree of recombination in the beam is likely small since there are no metal surfaces 
in the direct line-of-sight for the oxygen atoms as they enter the mass spectrometer 
ionizer. However, even if a large amount of recombination in the mass spectrometer does 
occur the mass spectrometer measurements would represent a lower limit on the 
percentage of atoms in the beam. So, the fast atom beam is conservatively estimated to 
be > 80 % atomic oxygen and < 20% molecular oxygen. Some oxygen ions are also 
present at less than 1 % of the beam. These ions are readily removed from the beam by 

the deflecnon plates described above. Impurities at other masses also represent less than 
1 to or the beam. 
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Summary Report on the Workshop on the Induced 
Environment of Space Station Freedom 

Huntsville, Alabama, January 30-31, 1991 

prepared by Dr Marsha R. Torr/MSFC 

Introduction, A workshop was held in Huntsville on January 30- 
SI, 1991 for the purpose of reviewing the state of knowledge of the 
likely induced external environment around Space Station Freedom. 
This workshop was chaired by Dr Marsha R. Torr and was a 
continuation of an activity coordinated by the Marshall Space Flight 
Center since 1987 and sponsored by the Office of Space Science and 
Applications and more recently by the Space Station Utilization 
Office at Reston. Two previous workshops have been held (one in 
October 1987 and one in May 1988) with the express purpose of 
assessing our understanding of the causative mechanisms underlying 
the various phenomena in the induced environment (glows, 
ionization, surface effects, gas envelope, etc). Both of the previous 
meetings led to the publication of NASA Conference Proceedings 
which document the contributions. As a result of the earlier 
reviews, a limited number of studies were funded in an attempt to 
obtain information on certain of the more fundamental processes 
involved. The activity has also been used to assist NASA 
Headquarters in assessing the impact of various Station design 
issues on potential attached payloads. Information gathered applies 
just as readily to critical Station systems such as the solar arrays. 

The purpose of this meeting was to address three questions: 

1) Where are we in our knowledge of the likely induced 
environment? 

2) What progress have we made in the past two years in 
understanding this environment? 

3) What areas of study are the most important for the next two 
years? 


Copies of the material presented at the meeting are available. 



Proceedings of Meeting. The agenda for the meeting is included 
here as Attachment I. The list of attendees is included as 
Attachment II. The meeting attendance was intentionally limited to 
the participants and NASA invitees. Some additional observers 
attended on a space available basis. The objective was to keep the 
meeting rather small so as to facilitate a good working environment. 
As with the previous two workshops, however, it was clear that as 
word of the meeting got out, the interest in attending was high and 
we had to turn a relatively large number of people away. This may 
indicate the desirability of scheduling a wider attendance "tutorial" 
meeting at some point to serve the purpose of briefing various 
personnel working on the Station design on the issues that this panel 
of experts is addressing. 

As part of the first question above, reviews were presented of 
various related activities that have been conducted over the past 
year, including a rather sizable DOD effort on vehicle glows, and a 
summary of the early findings from the LDEF program. In addition, 
presentations were made of the various studies that have been 
funded by this and other MSFC programs. 

Despite the very limited funding that has gone into the effort, the 
progress that has been made since the 1988 workshop was found to 
be considerable. This is largely due to the fact that the external 
Station environment represents a region of most unusual physical 
processes, including both gas phase and surface phase interactions, 
so that scientists in fields of appropriate expertize are challenged 
by the task of understanding the phenomena. A number of facilities 
have been put in place for the purpose of studying such environments 
in the laboratory and early results are beginning to emerge. The 
progress in developing models of the environment has made 
excellent progress. Almost three years ago we had very little to 
work with apart from unexplained Shuttle phenomena, and "back-of- 
the-envelope" projections for Space Station Freedom. Now we are at 
a point where we have the tools to begin to do believable 
evaluations. 



Recommendations. In order to predict the external Station 
environment, it is necessary that we are able to understand the 
nature of the processes taking place within it. At present, the state 
of knowledge is very poor. Over the past three years, with very 
nominal resources(~$150K per year), it has been possible to make 
progress in the field and to set in place important tools which are 
now available for use in further studies. 

At some point we anticipate that the Space Station Freedom will be 
used for purposes of attaching external payloads. If NASA is to 
remain a "smart buyer" in this area it is most important that we 
continue the efforts (albeit at a low level of effort) to understand 
the environment in which these systems will operate, and what the 
limitations on them will be. Maintenance of a multidiciplinary pool 
of expertise, such as is represented by the panels put together for 
these workshops and study efforts, is a valuable, if not essential, 
activity. 

At present there are several activities that we would recommend 
for immediate attention. For example, in order to conduct almost 
any assessment of the external environment, one must first model 
the gas concentrations surrounding the vehicle. Once this is known, 
one can proceed to model the induced optical emission and the 
induced ionization. With these established one can compute the 
backscatter fluxes, arcing, plasma discharges, optical thresholds, 
etc. One can optimize the placement of vents, and establish the 
impact of leaks. However, a fundamental input to models of the 
neutral gas environment is the gas phase collision cross sections. 
Those relevant to this problem are essentially unknown. It is 
important to measure these cross sections in the laboratory and 
then input them to the models. During a shuttle docking the Station 
surfaces can become flooded with contaminants. We do not know 
what the residence times of these contaminants are on the surfaces, 
and in the induced environment. These lifetimes can be measured in 
the laboratory and it is recommended that studies be made for key 
likely species such as water. If attached payloads are initially to be 
mounted on the JEM, it is important that we try to establish the 
environment for this location first. Relevant new evidence of 
synergistic effects is emerging from the LDEF studies. These need 
to be examined in the light of the Space Station Freedom external 
environment work. A summary of the specific recommendations 
made by individual panel members is included here as Attachment 
III. 



ATTACHMENT I: 


WORKSHOP on the INDUCED ENVIRONMENT of 
SPACE STATION FREEDOM 

Radisson Hotel, Huntsville, Alabama 

January 30-31, 1991 


8:30AM Introduction and Objectives of Workshop M. Torr 

8:45AM Environmental Definition and Assessment 

Program D. Brewer 

9:15AM Current Requirements and Plans to Verify 

These L. Leger 

9:45AM Status of Planned Investigations Attached to 

Space Station Freedom M. Sistilli 


10:15AM BREAK 


10:30AM Environmental Issues from the Utilization 

Viewpoint K. Schaefer 

11:00AM Vehicle Interaction Program at JHU/APL C. Meng 

1 1 :30 AM LUNCH 


1:00PM Glows, Accomodation and Surface Residence 

Times G. Caledonia 


1:30PM Summary of Space Station Grounding Issues R. Carruth 

2:00PM LDEF: Lessons Learned convenor: A. Whitaker 

R. Linton 
R. Rantanen 


3:15PM 


BREAK 



January 30. 1991. continued: 


3:30PM 

Overview of Spacecraft Glow 

T. Slanger 

4:00PM 

Summary of JPL Workshop on Modeling Tools 

J. Murphy 

4:30PM 

Discussion 


5:00PM 

End of Formal Presentations for the Day 


6:00PM 

SPECIAL EVENT AT ALABAMA SPACE AND ROCKET CENTER 
see attached sheet 

■Thursday. 

January 31. 1991 


8:30AM 

ISEM Space Station Model Update 

T. Gordon 

9:00 AM 

Cross Section Work at Rice 

K. Smith 

9:30AM 

Cross Section Work at UAH/MSFC 

C. Keffer 

10:00AM 

BREAK 


10:15AM 

Induced Emissions 

D. Torr 

10:45AM 

Atomic Oxygen Studies at PPPL J. Cuthbertson 

11:15AM 

Spacecraft Glow Studies and Mechanisms 

S. Mende 

1 1 :45AM 

LUNCH 


1 :00PM 

Visit to Space Station Mockup 


2:30PM 

Surface Temperature Dependent Glow 
Mechanisms in Space 

N. Tolk 

3:00PM 

What Have We Learned? 


3:30PM 

BREAK 


3:45PM 

4:30PM 

Where do we go from here? 
End of Meeting 




ATTACHMENT II 


LIST OF ATTTENDEES 


Marsha Torr, Chairperson 
Charles Keffer 
Anderson, Jeff 
Blanchette, Fred 
Brewer, Dana 
Caledonia, George 
Carruth, Ralph 
Chappell, Rick 
Collier, Jack 
Crane, Mike 
Cuthbertson, John 
Erlandson, Bob 
Espy, Pat 
Feddes, Allan 
Fichtl, George 
Gordon, Tim 
Hefling-Miller, Hilda 
Hwang, Kai-Shen 
Jongward, Gary 
Katz, Ira 
Leifer, Joel 
Linton, Roger 
McCombs , Roger 
Mende, Steve 
Melendez, Daniel 
Meng, Ching 
Murphy, Gerry 
Nahra, Henry 
Nebolsine, Peter 
O'Keefe, Ed 
Plaster, Teresa 
Rantanen, Ray 
Schaefer, Kevin 


NASA\MSFC 

UAH 

NASA/MSFC 

McDonnell Douglas Space Systems Co. 
NASA/HQ 

Physical Sciences Inc. 

NASA/MSFC 

NASA/MSFC 

NASA/HQ 

Boeing 

Princeton Plasma Physics Lab 

Applied Physics Lab 

NASA/MSFC 

Boeing 

NASA/MSFC 

Applied Science Technologies 

Grumman SSPSD 

Grumman SSEIC 

S- Cubed 

S-Cubed 

Booz-Allen & Hamilton 

NASA/MSFC 

BA&E 

Lockheed Palo Alto Research Lab 
NRC 

Applied Phyics Lab 
JPL 

NASA/ Lewis 

Physical Sciences Inc. 

Boeing 

Grumman Space Station Integration 

ROR Enterprises 

NASA/Reston 



Sistilli, Mark 

NASA /HQ 

SI anger, Tom 

SRI International 

Smith, Ken 

Rice University 

Snyder, Dave 

NASA /Lewis 

Suggs, Rob 

Grumman Space Station Integration 

Taylor, Bill 

NASA/HQ 

Tolk, Norman 

Vanderbilt University 

Torr, Doug 

UAH 

Whitaker, Ann 

NASA/MSFC 

Wyman, Pete 

NASA/ HQ 

Young, Dave 

Grumman 



ATTACHMENT III: 


Recommendations Made by Individual Panel Members (in alphabetical 
order). 


Brewer(HQ/MSS): 

‘Cross comparison of model results used by station engineers 
and science community using station vent, leakage and 
outgassing data for 

-each stage in assembly sequence to identify 
potential problems early on 

-identification of sensitive surfaces (which will 
change orientation as a function of stage) 

‘Identification/quantification of neutral/plasma effects on 
Station hardware. 

Carruth(MSFC/EH12) 

‘Needed: 

- definition of data required and models for environment 

- meeting of atomic oxygen facility groups 

- determine the causative processes involved in new 
LDEF phenomena (i.e. fluorescence and synergistic 
effects) 

- studies of basic interactions which affect surfaces 
rather subtly over time 

- investigations of plasma phenomena (interactions) on 
SSF systems and of mitigating techniques 

Cuthbertson(PPPL): 

For further progress we need: 

‘laboratory investigation of synergistic phenomena (UV, 
temperature cycling) 



Cross Section Work at UAH/MSFC 
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Radiation Damage Effects in Far Ultraviolet Filters and Substrates 

Charles E. Keffer 1 , Marsha R. Torr 2 , Muamer Zukic 1 , 

James F. Spann 2 , Douglas G. Torr 1 , and Jongmin Kim 1 

'Department of Physics 
The University of Alabama in Huntsville 
Optics Building, Suite 300 
Huntsville, AL 35899 

0 

2 N AS A/Marshall Space Flight Center 
Space Science Laboratory 
Huntsville, AL 35812 

ABSTRACT 

New advances in VUV thin film filter technology have been made using filter designs 
with multilayers of materials such as A1 2 0 3 , BaF 2 , CaF 2 , FFfO z , LaF 3 , MgF 2 , and Si0 2 . Our 
immediate application for these filters will be in an imaging system to be flown on a satellite 
where a 2 X 9 R E orbit will expose the instrument to approximately 275 krads of radiation. In 
view of the fact that no previous studies have been made on potential radiation damage of these 
materials in the thin film format, we report on such an assessment here. Transmittances and 
reflectances of BaF 2 , CaF 2 , Hf0 2 , LaF 3 , MgF 2 , and Si0 2 thin films on MgF 2 substrates, A1 2 0 3 thin 
films on fused silica substrates, uncoated fused silica and MgF 2 , and four multilayer filters made 
from these materials were measured from 120 nm to 180 nm before and after irradiation by 250 
krads from a “Co gamma radiation source. No radiation-induced losses in transmittance or 
reflectance occurred in this wavelength range. Additional postradiation measurements from 160 
nm to 300 nm indicated a 3 - 5% radiation- induced absorption near 260 nm in some of the 
samples with MgF 2 substrates. From these measurements it is concluded that far ultraviolet 
filters made from the materials tested should experience less that 5 % change from exposure to 
up to 250 krads of high energy radiation in space applications. 

1. INTRODUCTION 

The Ultraviolet Imager (UVI) for the Global Geospace Science (GGS) POLAR spacecraft 
has been designed to acquire coherent global images of the Earth’s aurora under both day and 
night conditions 1 . The major goals of the observations are to obtain information on the temporal 
and spatial morphology of the aurora, as well as on the total incident energy flux and energy 
characteristics of the precipitating particles giving rise to the aurora. In order to achieve these 
goals, a fast (f/2.9) imaging camera has been developed which obtains wide field of view (8°) 
images through 5 specially designed FUV filters. Success in achieving the scientific goals is 
therefore dependent on the high performance of these filters. In addition to spectral resolution 
of closely spaced features (130.4 nm and 135.6 nm), the filters must also have sufficiently high 
throughput in order to allow detection of the relatively weak Lyman Birge Hopfield bands of N 2 , 
and must be capable of blocking out-of-band sunlight by a factor of 10 4 at wavelengths longer 
than 190 nm and isolating the OI emissions at 130.4 nm and 135.6 nm from the bright Lyman-a 
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emission at 121.6 nm. Attainment of these filter characteristics has required the development of 
a new mu n ayer stack named the II multilayer. 2 The design approach was to obtain high filter 
throughput by minimizing energy losses in the filters with the II multilayer design and bv 
operating in a reflective, rather than transmissive, mode. Each UVI filter consists of three II 
multilayer reflection filters and one transmission filter in series to achieve the necessarv 
performance requirements. The transmission filter provides short wavelength rejection by acting 

“ * , CUt n, I!" Wh / e ? e reflectors provide excellent long wavelength rejection 8 
At ^ ?Tvi ^ 0Cking °£ s 1 unll .f ht beyond 190 nm is provided by a solar blind Csl photocathode 
in the UVI detector. A detailed description of the design and development of the UVI filters is 
given by Zukic, et al. 

The three year mission lifetime of the POLAR spacecraft will result in a radiation 
exposure to the UVI of approximately 275 krads. Radiation damage to optical materials has a 
long history of study m the space science and nuclear physics fields. At ultraviolet wavelengths 
transmission of optical materials exposed to high energy radiation has been previously measured 
SWISS' ° fbuUc materials including A1 2 0«, BaFj 4 * 7 , CaF 2 4 , CsF 6 , fused silica 4 - 5 - 8 * 9 LiF 4 - 510 -’ 1 
MgF 2 ■ and several optical glass materials 4 - 8 . Ultraviolet transmission has’ also been 
measured for irradiated colored glass filters and thin film interference filters 9 . Ultraviolet 
reflectance measurements of irradiated optical materials have been made at 121.6 nm for 
aluminum mirrors with MgF 2 protective coatings 511 . Radiation damage effects on the materials 
used m these FUV filters have not been previously studied in the thin film format. In this paper 
t 3 n results of 111 investigation into the potential for radiation induced damage to the 
UVI filters by making measurements on the thin films and substrates used to fabricate the filters. 
Only long term radiation damage was studied since measurements were not made immediately 
following radiation exposure. 


2. EXPERIMENTAL PROCEDURES 

A set of 22 samples was selected for use in this study. These include BaF 2 , CaF HfO 
LaF 3 , MgF 2 , and Si0 2 single layer thin films on MgF 2 substrates (two of each), A1 2 0 3 single layer 
thm films on fused silica substrates (two), two uncoated MgF 2 substrates, and two uncoated 
fused silica substrates. Two multilayer transmission filters (BaFj/MgF 2 on a MgF 2 substrate) and 
two multilayer reflection filters (LaFj/MgFj on a fused silica substrate) were also used Every 
single layer thin film has a thickness between 8 nm and 22 nm with the exact thickness for each 
film material corresponding to the optimum thickness from the II multilayer design. The filters 
used are 35-layer II stacks. The complete set of samples is summarized in Table 1 and is 
representative of the substrates, thin films, and filters which make up the UVI filter system. 

All of the test samples were prepared with 12.7 mm diameter and 2 mm thick VUV grade 
substrates used as received from Acton Research Corporation. BaF 3 , CaF 3 and LaF film 
materials were 99.9% pure from CERAC. A1 2 0 3 (99.5%), MgF 2 (99.95%), and Si0 2 (99?98%) 
thin films were deposited with standard Balzers coating materials. Hf0 2 was supplied by EM 
Chemicals with 99.5% purity. All depositions were completed at The University of Alabama in 
Huntsville Optical Aeronomy Laboratory. An oil-free coating chamber with a sorption pump and 
cryopump was used to prevent hydrocarbon contamination of the films. 


SPIE Vol. 2018 Passive Materials for Optical Elements II (1993) / 93 



f ,u Ol I e . 0f eac i! !ype of Slngle layer 111111 fihn ’ one of each type Of uncoated substrate and all 
o e multilayer filters were mounted in an aluminum sample holder which was then sealed in 

an ultrahigh punty dry nitrogen purged double bag made from ultra clean radiation resistant 

/T WCre S “ bsequentI y ^ to Goddard s Pace Right Center where they were 

exposed for 1 10 hours to a Co gamma radiation source. The total dose received was 250 krads 
After radiation exposure, the Aclar bag remained clear and had no visible evidence of radiation 
damage. The remaining single layer thin films and uncoated substrates served as control samples 
and received identical handling, with the exception of radiation exposure, including being sealed 
in an ultrahigh punty dry nitrogen purged double bag made from Aclar. 


c ?? ° f rcflectance and transmittance measurements was performed at the Marshall 

bpace Right Center Ionosphere Thermosphere Mesosphere Branch on the 22 samples described 
above. The preradiation measurements were performed in a hydrocarbon-free cryopumped 
vacuum chamber at a pressure below 10 5 Torr. A deuterium lamp with a MgF 2 window 

Wlth vacuum monochromator provided 1 nm FWHM spectral resolution over 

the 120 nm to 180 nm wavelength range. Folding and collimating optics and a 6 mm diameter 
aperture limited the light incident on the eight position filter wheel holding the 12.7 mm diameter 
substrates to an area approximately l/4th the area of the thin film. The small aperture insured 
that small positioning errors from rotating the filter wheel did not significantly impact the 
reproducibility of the measurements. A photomultiplier tube with a MgF 2 window and a semi- 
transparent Cs-I photocathode served as the detector for all measurements. The postradiation 
measurements were accomplished in a different hydrocarbon-free cryopumped vacuum chamber 
also at a pressure below 10 Torr. A similar deuterium lamp, 0.2-m vacuum monochromator, 
and aperture were used for the postradiation measurements. However, the second chamber did 
not have any folding and collimating optics. The same Cs-I photomultiplier tube and eieht 
position filter wheel was used for both the pre- and postradiation measurements in the 120 nm 
to 180 nm wavelength range. A second set of postradiation transmission measurements was 
made from 160 nm to 300 nm using a Cs-Te photomultiplier tube. All reflectance measurements 
were performed at a 45° angle of incidence while transmittance was measured at normal 
incidence. Figures 1 and 2 illustrate the optical configurations for the two vacuum chambers 


For the reflection measurements, two identical scans from 120 nm to 180 nm were 
performed on every sample both before and after the radiation exposure. A reference reflection 
fiiterwas measured as a control with every reflection measurement to check the reliability and 
reproducibility of the data from measurement to measurement. The reference filter also insured 
that no systematic errors were introduced when the vacuum chamber was vented and a new set 
of samples was installed. Two identical transmission measurements from 120 nm to 180 nm 
were also performed on every sample both before and after the radiation exposure. An uncoated 
fused silica substrate in the incident beam served as a filter to block second order monochromator 
reflections for the 160 nm to 300 nm postradiation transmission measurements. Two identical 
scans from 160 nm to 300 nm were recorded for each sample. A reference transmission filter 
was measured with every transmission measurement as a check on the measurement reliability 
and reproducibility. The reflection and transmission reference filters also provided traceability 
between the two vacuum chambers used in the preradiation and postradiation measurements. 
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fmm Ju xt ^ V ^ UCS were de£ermined b y comparing the reflected beam intensity 

from each filter with the reflected beam intensity from a VUV-enhanced aluminum mirror whose 

ectance was calibrated using an Acton vacuum reflectometer model VRSC-100 Absolute 

rrx« Val , U H S K W ' re “““ by meaSUri " S rad0 of “am“ 

r: Md si * nai drift rr 


3. RESULTS 

The primary wavelength range of interest for the UVI filters is 120 nm to 180 nm 
Transmission and reflection measurements in this range were analyzed by first averaging the two 
for each sample. The difference, defined as the average of the 
measurements nunus the average of the preradiarion measurements, was then calculated. Figures 

“ U . gh , U . show results for a]1 of films studied. The error bars represent one 

f eVia “°"- The prcradianon measurements have a lower signal to noise ratio due to 

J" intensi£ y losses from the folding and collimating optics used in the preradiation 
measurements but not in the postradiation measurements. 

Repeated measurements of the transmission and reflection reference filters result in an 
esnmated measurement uncertainty (a) of 2%. Differences between the postradiation and 

lZr^Tr t m T UTem ^ tS ° f * rea “ r *“ 4% 316 considercd to be due to radiation-induced 
e tects. Particular regard is given to those wavelengths where radiation-induced absorption bands 

are known to occur in the bulk materials. 4 With this criterion as a benchmark, none of the 13 
samples which each received 250 krads of radiation showed evidence of statistically significant 
radianon-induced changes in reflectance in the 120 nm to 180 nm wavelength range For the 

" leaSUremen£s in the 120 nm t0 180 "nr range, only the MgF 2 uncoated substrate 
had any statistically significant changes. This sample had an apparent absorption short of 160 
nm before radiation exposure which was not present after 250 loads radiation exposure. 

Most of the UV optical materials that have been studied in bulk form have their strongest 
radianon-induced absorpnon bands at wavelengths longer than 180 nm. 4 From 160 nm to 300 
nm, die irradiated samples were compared with their equivalent control samples since no 
preradianon measurements were made in this wavelength range. Differences were calculated in 
the same manner as for the 120 nm to 180 nm wavelength measurements. Of the nine irradiated 

* n thlS S UP ’ th ^ C Sh ° W cvldence of a sma11 radiation-induced absorption band centered 
near 250 nm -260 nm. Figures 5, 8, and 9 are for a CaF 2 film on a MgF 2 substrate, a LaF 3 film 
on a MgF 2 substrate, and a MgF 2 film on a MgF 2 substrate, respectively. The coincidence of the 
absorption with the well known 260 nm absorption band in MgF 2 4 and the fact that the MgF, 

c “ t0 thC sam P les strc >ngly suggests that the radiation damage is limited 
Ae ^ g f 2 buUc matenal - The small 3 - 5% absorption seen in these three samples is near the 
limit of the measurement uncertainty and thus explains how four other samples with MgF, 
substrates do not show similar indications of radiation-induced absorption. 

M,v f n ea,h /T d toadiated a >' 5 thiak MgF 2 aystal with 10” electrons/cm 1 at I 

McV followed by 10 electrons/cm 1 at 2 MeV. They observed 16% loss in transmittance at 120 
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nm and approximately 60% loss in transmittance at 260 nm. Hass and Hunter 5 reported a 5% 
decrease in transmittance at 121.6 nm and approximately 20% decrease at 260 nm after exposure 
of their 3 mm thick MgF 2 crystal to 10 15 electrons/cm 2 at 1 MeV. They attribute the smaller loss 
m transmittance compared with Heath and Sacher’s results to a MgF 2 crystal which was probably 
higher in punty. Becher, et al. found approximately 20% decrease in MgF 2 transmittance at 260 
nm from ireadianon with 8 X 10 protons/cm 2 at 85 MeV. Reft, et al. u measured MgF 
transmittance losses of 10% at 121.6 nm and 14% at 180 nm from irradiation with 85 - and 600 2 
MeV protons at a total absorbed energy of 5.7 X 10 14 MeV/cm 3 . In the present study, the total 
absorbed energy in the samples with MgF 2 substrates is 5 X 10 13 MeV/cm 3 . This is over ten 
times less total energy absorbed than in any of the other reports cited. The maximum 5% 
transmittance loss at 260 nm and no measurable change at 120 nm reported in the present study 

is thus in reasonable agreement with the earlier reports on radiation-induced damage in MgF 2 
crystals. 


Fused silica has been reported by Heath and Sacher 4 to decrease in transmittance 
approximately 15% at 220 nm after irradiation with 10 14 electrons/cm 2 at 1 MeV and 10 14 
ele ~™ ns/cm r at 2 MeV - Hass and Hunter 5 measured approximately 10% loss in transmittance 
at 220 nm after irradiation of fused silica by 10 15 electrons/cm 2 at 1 MeV. Nicoletta and 
ubanks used 10 electrons/cm 2 at 1.5 MeV to irradiate fused silica and reported similar 
transmittance losses to Heath and Sacher. In the present study, no radiation-induced transmission 
losses were observed in the samples with fused silica substrates. The total absorbed energy was 
3 X 10' 3 MeV/cm 3 for these samples. 


, . Nicoletta and Eubanks measured radiation damage in three interference filters which were 

shielded by 3. i mm of fused silica and then inadiated with 2.7 X 10 13 electrons/cm 2 of energies 

0. 3 MeV, 0.5 MeV, 1.0 MeV, and 1.5 MeV. Transmission losses were very small from the 
electron irradiation shielded by the fused silica. No changes in transmittance or reflectance have 
been observed in the present study which can be attributed to radiation damage in the thin films 
used. This is probably because only a small fraction of the energy is likely to be absorbed in 
the thin films and also because the radiation levels used were moderate. 

4. CONCLUSION 

None of the samples tested showed evidence of any major radiation-induced changes in 
reflectance or transmittance. Since the UVI filters are made from m ate rials which we have tested 
and smee they are housed within the body of the instrument where they will experience a total 
radiation dose of less that 20 krads, it is anticipated that they will suffer negligible loss in 
reflectance or transmittance during the nominal 3 year mission lifetime of the POLAR spacecraft. 
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Table 1: Summary of thin films and substrates tested. 


Film Type 

Film Materials 

Substrate 

Single Layer 

A1A 

Fused Silica 

Single Layer 

BaF 2 

MgF 2 

Single Layer 

CaF 2 

MgF 2 

Single Layer 

Hf0 2 

MgF 2 

Single Layer 

LaF 3 

MgF 2 

Single Layer 

MgF 2 

MgF 2 

Single Layer 

Si0 2 

MgF 2 

n Multilayer 

BaF 2 /MgF 2 

MgF 2 

11 Multilayer 

LaFj/MgF 2 

Fused Silica 

None 

None 

Fused Silica 

None 

None 

MgF 2 



Figure 1: Vacuum chamber for preradiation measurements. Chamber is hydrocarbon free and 
cryogenically pumped. Optical configuration shown is for 45° reflectance measurements. 
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Figure 2: Vacuum chamber for postradiation measurements. Chamber is hydrocarbon free and 
cryogenicaUy pumped. Optical configuration shown is for 45° reflectance measurements using 
a Cs-I pmt for the 120 nm to 180 nm range. 
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Figure 3: Reflectance and transmittance of an A1 2 0 3 thin film on a fused silica substrate. 
Comparison is shown in the 120 nm to 180 nm range before and after radiation exposure of a 
single sample. In the 160 nm to 300 nm range, comparison is between the radiation exposed 
sample and a control sample with no radiation exposure. 
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Figure 4: Reflectance and transmittance of a BaF 2 thin film on a MgF 2 substrate. Comparison 
is shown in the 120 nm to 180 nm range before and after radiation exposure of a single sample. 
In the 160 nm to 300 nm range, comparison is between the radiation exposed sample and a 
control sample with no radiation exposure. 



120 140 160 ISO 200 220 240 260 280 300 


Wavelength (nm) 

Figure 5: Reflectance and transmittance of a CaF 2 thin film on a MgF 2 substrate. Comparison 
is shown in the 120 nm to 180 nm range before and after radiation exposure of a single sample. 
In the 160 nm to 300 nm range, comparison is between the radiation exposed sample and a 
control sample with no radiation exposure. 
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Figure 6: Reflectance and transmittance of an uncoated fused silica substrate. Comparison is 
shown in the 120 nm to 180 nm range before and after radiation exposure of a single sample, 
in the 160 nm to 300 nm range, comparison is between the radiation exposed sample and a 
control sample with no radiation exposure. 
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Figure 7: Reflectance and transmittance of a Hf0 2 thin film on a MgF 2 substrate. Comparison 
is shown in the 120 nm to 180 nm range before and after radiation exposure of a single sample. 
In the 160 nm to 300 nm range, comparison is between the radiation exposed sample and a 
control sample with no radiation exposure. 
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Figure 8: Reflectance and transmittance of a LaF 3 thin film on a MgF 2 substrate. Comparison 

In rh 0 T™ thC 12 °JT t0 180 nm 131166 before and ^ radiation exposure of a single sample. 
In the 160 nm to 300 nm range, comparison is between the radiation exposed sample and a 
control sample with no radiation exposure. 
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igure 9: Reflectance and transmittance of a MgF 2 thin film on a MgF 2 substrate. Comparison 
is shown in the 120 nm to 180 nm range before and after radiation exposure of a single sample, 
in the 160 nm to 300 nm range, comparison is between the radiation exposed sample and a 
control sample with no radiation exposure. 
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150 nm to 300 nm range, comparison is between the radiation exposed sample and a control 
sample with no radiation exposure. 
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Figure 1 1: Reflectance and transmittance of a Si0 2 thin film on a MgF 2 substrate. Comparison 
is s hown in the 120 Iran to 180 nm range before and after radiation exposure of a single sample. 

m the 160 nm to 300 nm range, comparison is between the radiation exposed sample and a 
control sample with no radiation exposure. 
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Abstract 


New advances in VUV thin film filter technology have been made using filter designs 
with multilayers of materials such as A1A, BaF 2 , CaF 2 , Hf0 2 , LaF 3 , MgF 2 , and Si0 2 . Our 
immediate application for these filters will be in an imaging system to be flown on a satellite 
where a 2 X 9 R e orbit will expose the instrument to approximately 250 krads of radiation. In 
view of the fact that no previous studies have been made on potential radiation damage of these 
materials in the thin film format, we report on such an assessment here. Transmittances and 
reflectances of BaF 2 , CaF 2 , Hf0 2 , LaF 3 , MgF 2 , and Si0 2 thin films on MgF 2 substrates, A1 2 0 3 thin 
films on fused silica substrates, uncoated fused silica and MgF 2 , and four multilayer filters made 
from these materials were measured from 120 nm to 180 nm before and after irradiation by 250 
krads from a “Co gamma radiation source. No radiation-induced losses in transmittance or 
reflectance occurred in this wavelength range. Additional postradiation measurements from 160 
nm to 300 nm indicated a 2 - 5% radiation-induced absorption near 260 nm in some of the 
samples with MgF 2 substrates. From these measurements it is concluded that far ultraviolet 
filters made from the materials tested should experience less that 5% change from exposure to 
up to 250 krads of high energy radiation in space applications. 

Key words : Radiation effects, ultraviolet filters, thin films 

A. Introduction 

The Ultraviolet Imager (UVI) for the Global Geospace Science (GGS) POLAR spacecraft 
has been designed to acquire coherent global images of the Earth’s aurora under both day and 
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night conditions*. The major goals of the observations are to obtain information on the temporal 
and spatial morphology of the aurora, as well as on the total incident energy flux and energy 
characteristics of the precipitating particles giving rise to the aurora. In order to achieve these 
goals, a fast (f/2.9) imaging camera has been developed which obtains wide field of view (8°) 
images through 5 specially designed FUV filters. Success in achieving the scientific goals is 
therefore dependent on the high performance of these filters. In addition to spectral resolution 
of closely spaced features (130.4 nm and 135.6 nm), the filters must also have sufficiently high 
throughput in order to allow detection of the relatively weak Lyman Birge Hopfield bands of N 2 , 
and must be capable of blocking out-of-band sunlight by a factor of 10* at wavelengths longer 
than 190 nm and isolating the OI emissions at 130.4 nm and 135.6 nm from the bright Lyman-a 
emission at 121.6 nm. Attainment of these filter characteristics has required the development of 
a new multilayer stack named the n multilayer. 2 The design approach was to obtain high filter 
throughput by minimizing energy losses in the filters with the II multilayer design and by 
operating in a reflective, rather than transmissive, mode. Each UVI filter consists of three II 
multilayer reflection filters and one transmission filter in series to achieve the necessary 
performance requirements. The transmission filter provides short wavelength rejection by acting 
as a cut-on filter while the three reflectors provide excellent long wavelength rejection. 
Additional blocking of sunlight beyond 190 nm is provided by a solar blind Csl photocathode 
in the UVI detector. A detailed description of the design and development of the UVI filters is 
given by Zukic, et al. 3 

The three year mission lifetime of the POLAR spacecraft will result in a radiation 
exposure to the UVI of approximately 250 loads due to passage of the spacecraft through the 
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Earth’s radiation belts. This radiation environment consists primarily of electrons and protons 
with energies from less than 1 MeV to tens of MeV. Bremsstrahlung from high energy electrons 
impacting the spacecraft or instrument housings can also be a major source of radiation. We 
have taken the common approach of using a ®°Co gamma radiation source to simulate the Earth’s 
radiation environment In order for this to be a valid approach for radiation damage assessments 
in optical materials, the density of F-center formation from different types of radiadon must be 
known. Durand, et al. 4 have found that the density of F-center formation from x-ray 
bombardment of LiF is proportional to the square root of the total energy absorbed. Becher, et 
al. 5 have found that the same relationship holds true for proton and electron bombardment of LiF 
and MgF 2 with some deviation below a few Mrad in MgF 2 . The use of a "Co radiation source 
to simulate radiation damage to optical materials in space should be a reasonable approximation 

particularly for optical elements which are fairly will shielded form direct exposure to the proton 
and electron flux. 

Radiation damage to optical materials has a long history of study in the space science and 
nuclear physics fields. At ultraviolet wavelengths, transmission of optical materials exposed to 
high energy radiation has been previously measured for a variety of bulk materials including 
A1 2 Oj , BaFj 6-9 , CaF 2 6 , CsF*, fused silica 6 ’ 7,10 ’", LiF 5 ' 7,12 , MgF^ 7,12 , and several optical glass 
materials 6,10 . Ultraviolet transmission has also been measured for irradiated colored glass filters 
and thin film interference filters". Ultraviolet reflectance measurements of irradiated optical 
materials have been made at 121.6 nm for aluminum mirrors with MgF 2 protective coatings 7,12 . 
Radiation damage effects on the materials used in the UVI filters have not been previously 
studied in the thin film format. In this paper we report the results of an investigation into the 
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potential for radiation-induced damage to the UVI filters by making measurements on the thin 
films and substrates used to fabricate the filters. Only long term radiation damage was studied 
since measurements were not made immediately following radiation exposure. 

B. Experimental Procedures 

A set of 22 samples was selected for use in this study. These include BaF 2 , CaF 2 , Hf0 2 , 
LaF 3 , MgF 2 , and Si0 2 single layer thin films on MgF 2 substrates (two of each), A1 2 0 3 single layer 
thin films on fused silica substrates (two), two uncoated MgF 2 substrates, and two uncoated 
fused silica substrates. Two multilayer transmission filters (BaFj/MgF 2 on a MgF 2 substrate) and 
two multilayer reflection filters (LaFj/MgF 2 on a fused silica substrate) were also used. Every 
single layer thin film has a thickness between 8 nm and 22 nm with the exact thickness for each 
film material corresponding to the optimum thickness from the n multilayer design. The filters 
used are 35-layer II stacks. The complete set of samples is representative of the substrates, thin 
films, and filters which make up the UVI filter system. 

All of the test samples were prepared with 12.7 mm diameter and 2 mm thick VUV grade 
substrates used as received from Acton Research Corporation. BaF 3 , CaF 3 , and LaF 3 film 
materials were 99.9% pure from CERAC. A1 2 0 3 (99.5%), MgF 2 (99.95%), and SiOj (99.98%) 
thin films were deposited with standard Balzers coating materials. Hf0 2 was supplied by EM 
Chemicals with 99.5% purity. All depositions were completed at The University of Alabama in 
Huntsville Optical Aeronomy Laboratory. An oil-free coating chamber with a sorption pump and 
cryopump was used to prevent hydrocarbon contamination of the films. 

One of each type of single layer thin film, one of each type of uncoated substrate and all 


5 



of the multilayer filters were mounted in an aluminum sample holder which was then sealed in 
an ultrahigh purity dry nitrogen purged double bag made from ultra clean radiation resistant 
Aclar. The samples were subsequently taken to Goddard Space Flight Center where they were 
exposed for 1 10 hours to a “Co gamma radiation source. The total dose received was 250 krads. 
After radiation exposure, the Aclar bag remained clear and had no visible evidence of radiation 
damage. The remaining single layer thin films and uncoated substrates served as control samples 
and received identical handling, with the exception of radiation exposure, including being sealed 
in an ultrahigh purity dry nitrogen purged double bag marir. from Aclar. 

A series of reflectance and transmittance measurements was performed at the Marshall 
Space Flight Center Ionosphere Thermosphere Mesosphere Branch on the 22 samples described 
above. The preradiation measurements were performed in a hydrocarbon-free cryopumped 
vacuum chamber at a pressure below 1CT 5 Torr. A deuterium lamp with a MgF 2 window 
together with a 0.2-m vacuum monochromator provided 1 run FWHM spectral resolution over 
the 120 nm to 180 nm wavelength range. Folding and co llimatin g optics and a 6 mm diameter 
aperture limited the light incident on the eight position filter wheel holding the 12.7 mm diameter 
substrates to an area approximately l/4th the area of the thin film. The small aperture insured 
that small positioning errors from rotating the filter wheel did not significantly impact the 
reproducibility of the measurements. A photomultiplier tube with a MgF 3 window and a semi- 
transparent Cs-I photocathode served as the detector for all measurements. The postradiation 
measurements were accomplished in a different hydrocarbon-free cryopumped vacuum chamber 
also at a pressure below 10 5 Torr. A similar deuterium lamp, 0.2-m vacuum monochromator, 
and aperture were used for the postradiation measurements. However, the second chamber did 
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not have any folding and collimating optics. The same Cs-I photomultiplier tube and eight 
position filter wheel was used for both the pre- and postradiation measurements in the 120 nm 
to 180 nm wavelength range. A second set of postradiation transmission measurements was 
made from 160 nm to 300 nm using a Cs-Te photomultiplier tube. All reflectance measurements 

were performed at a 45° angle of incidence while transmittance was measured at normal 
incidence. 

For the reflection measurements, two identical scans from 120 nm to 180 nm were 
performed on every sample both before and after the radiation exposure. A reference reflection 
filter was measured as a control with every reflection measurement to check the reliability and 
reproducibility of the data from measurement to measurement. The reference filter also insured 
that no systematic errors were introduced when the vacuum chamber was vented and a new set 
of samples was installed. Two identical transmission measurements from 120 nm to 180 nm 
were also performed on every sample both before and after the radiation exposure. An uncoated 
fused silica substrate in the incident beam served as a filter to block second order monochromator 
reflections for the 160 nm to 300 nm postradiation transmission measurements. Two identical 
scans from 160 nm to 300 nm were recorded for each sample. A reference transmission filter 
was measured with every transmission measurement as a check on the measurement reliability 
and reproducibility. The reflection and transmission reference filters also provided traceability 
between the two vacuum chambers used in the preradiation and postradiation measurements. 

Absolute reflectance values were determined by comparing the reflected beam intensity 
from each filter with the reflected beam intensity from a VUV -enhanced al uminum minor whose 
reflectance was calibrated using an Acton vacuum reflectometer model VRSC-100. Absolute 
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transmittance values were determined by measuring the ratio of the transmitted beam intensity 
to the unattenuated beam intensity measured at an empty filter wheel position. Background 
signal was subtracted and signal drift was corrected using measurements before and after the 
filters were measured. 

C. Results 

The primary wavelength range of interest for the UVI filters is 120 nm to 180 am. 
Transmission and reflection measurements in this range were analyzed by averaging the two 
scans taken for each sample. Repeated measurements of the transmission and reflection reference 
filters result in an estimated measurement uncertainty (o) of ± 2%. Differences between the 
postradiation and preradiation measurements of greater than 2% are considered to be due to 
radiation-induced effects. Particular regard is given to those wavelengths where radiation-induced 
absorption bands are known to occur in the bulk materials. 6 With this criterion as a benchmark, 
none of the 13 samples which each received 250 krads of radiation showed evidence of 
statistically significant radiation-induced changes in reflectance in the 120 nm to 180 nm 
wavelength range. For the transmittance measurements in the 120 nm to 180 nm range, only the 
MgF 2 uncoated substrate had any statistically significant changes. This sample had an apparent 
absorption short of 160 nm before radiation exposure which was not present after 250 krads 
radiation exposure (see Figure 1). 

Most of the UV optical materials that have been studied in bulk form have their strongest 
radiation- induced absorption bands at wavelengths longer than 180 nm. 6 From 160 nm to 300 
nm, the irradiated samples were compared with their equivalent control samples since no 
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preradiation measurements were made in this wavelength range. Of the nine irradiated samples 
in this group, only the CaF 2 film on a MgF 2 substrate (Figure 2) and the LaF 3 film on a MgF 2 
substrate (Figure 3) show evidence of a small radiation-induced absorption band centered near 
250 nm - 260 nm. The coincidence of the absorption with the well known 260 nm F-center 
absorption band in MgF 2 6 and the fact that the MgF 2 substrate is common to these samples 
strongly suggests that the radiation damage is limited to the MgF 2 bulk material. The small 2 - 
5% absorption seen in these samples is near the limit of the measurement uncertainty and thus 
explains how the other samples with MgF 2 substrates do not show similar indications of 
radiation-induced absorption. 

The four multilayer filter samples were measured in the wavelength range from 120 nm 
to 180 nm. None of these filters show any statistically significant radiation-induced changes in 
transmittance or reflectance. Table 1 summarizes the results for all the filters, thin films, and 
substrates tested. 

Heath and Sacher 6 irradiated a 1.5 mm thick MgF 2 crystal with 10 14 electrons/cm 2 at 1 
MeV followed by 10 14 electrons/cm 2 at 2 MeV. They observed 16% loss in transmittance at 120 
nm and approximately 60% loss in transmittance at 260 nm. Hass and Hunter 7 reported a 5% 
decrease in transmittance at 121.6 nm and approximately 20% decrease at 260 nm after exposure 
of their 3 mm thick MgF 2 crystal to 10 15 electrons/cm 2 at 1 MeV. They attribute the smaller loss 
in transmittance compared with Heath and Sacher’ s results to a MgF 2 crystal which was probably 
higher in purity. Becher, et al. s found approximately 20% decrease in MgF 2 transmittance at 260 
nm from irradiation with 8 X 10 12 protons/cm 2 at 85 MeV. Refit, et al. 12 measured MgF 2 
transmittance losses of 10% at 121.6 nm and 14% at 180 nm from irradiation with 85 - and 600- 
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MeV protons at a total absorbed energy of 5.7 X 10 14 MeV/cm 3 . In the present study, the total 
absorbed energy in the samples with MgF 2 substrates is 5 X 10 13 MeV/cm 3 . This is over ten 
times less total energy absorbed than in any of the other reports cited. The mavimum 5% 
transmittance loss at 260 nm and no measurable change at 120 nm reported in the present study 
is thus in reasonable agreement with the earlier reports on radiation-induced damage in MgF 2 
crystals. 

Fused silica has been reported by Heath and Sacher 6 to decrease in transmittance 
approximately 15% at 220 nm after irradiation with 10 14 electrons/cm 2 at 1 MeV and 10 14 
electrons/cm 2 at 2 MeV. Hass and Hunter 7 measured approximately 10% loss in transmittance 
at 220 nm after irradiation of fused silica by 10 15 electrons/cm 2 at 1 MeV. Nicoletta and 
Eubanks" used 10 14 electrons/cm 2 at 1.5 MeV to irradiate fused silica and reported similar 
transmittance losses to Heath and Sacher. 6 In the present study, no radiation-induced 
transmission losses were observed in the samples with fused silica substrates. The total absorbed 
energy was 3 X 10 13 MeV/cm 3 for these samples. 

Nicoletta and Eubanks" measured radiation damage in three interference filters which 
were shielded by 3.1 mm of fused silica and then irradiated with 2.7 X 10 13 electrons/cm 2 of 
energies 0.3 MeV, 0.5 MeV, 1.0 MeV, and 1.5 MeV. Transmission losses were very small from 
the electron irradiation shielded by the fused silica. No changes in transmittance or reflectance 
have been observed in the present study which can be attributed to radiation damage in the thin 
films used. This is probably because only a small fraction of the energy is likely to be absorbed 
in the thin films and also because the radiation levels used were moderate. 
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D. Conclusion 


None of the samples tested showed evidence of any major radiation- induced changes in 
reflectance or transmittance. Since the UVT filters are made from materials which we have tested 
and since they are housed within the body of the instrument where they will experience a total 
radiation dose of less that 20 loads, it is anticipated that they will suffer negligible loss in 
reflectance or transmittance during the nominal 3 year mission lifetime of the POLAR spacecraft. 
Thus, we have experimentally demonstrated that II multilayer thin films are sufficiently radiation 
resistant to allow their use in the high energy radiation environment of space. In a follow-up 
study, we will present a theoretical model of high energy radiation effects on multilayer thin 
films. 

This work was supported by NASA grant NAG8-834 and NASA contract NAS8-38145. 
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Figure Captions 


Figure 1: Transmission measurements for the uncoated MgF 2 substrate. Comparison is shown 
in the 120 nm to 180 nm range before and after 250 krads of radiation exposure to the sample. 
The error bars represent one standard deviation. 

Figure 2: Transmission measurements for the CaF 2 thin film on a MgF 2 substrate. Comparison 
in the 160 nm to 300 nm range is between the radiation exposed sample and a control sample 
with no radiation exposure. The error bars represent one standard deviation. 

Figure 3: Transmission measurements for the LaF 3 thin film on a MgF 2 substrate. Comparison 
in the 160 nm to 300 nm range is between the radiation exposed sample and a control sample 
with no radiation exposure. The error bars represent one standard deviation. 
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Use of FUV auroral emissions as diagnostic indicators 

G. A. Germany , 1 M. R. Torr , 2 D. G. Torr , 3 and P. G. Richards 4 

In an earlier study we modeled selected FUV auroral emissions (O I (1356 A), N 2 Lyman -Birge-Hopfield 
(LBH) (1464 A), and LBH (1838 A)) to examine the sensitivity of these emissions and their ratios to likely 
changes in the neutral atmosphere. In this paper we extend that study to examine the dependence of these 
same emissions and their ratios on the shape of the energy distribution of the auroral electrons. In particular, 
we wish to determine whether changes in energy spectra might interfere with our determination of the 
characteristic energy. Modeled column -integrated emissions show relatively small (<30%) dependences on 
the shape and width of the incident energy spectrum, provided the average energy and total energy flux of 
the energy distribution are held constant. Long -wavelength FUV emissions, which are relatively unaffected 
by 0 2 absorption losses, exhibit virtually no dependence on the shape of the incident energy distribution. 
Changes in ratios of FUV short- to long-wavelength emissions as a function of characteristic energy are 
much larger than those due to changes in energy distribution. As a result, the determination of characteristic 
energy using these emission ratios is relatively unambiguous. We also examine the relative intensities of the 
aurora and the dayglow for various conditions. The intensities of modeled FUV auroral emissions relative to 
the dayglow emissions are presented as a function of solar zenith angle and incident energy flux. Under 
certain conditions (energy flux £ 1 erg enr 2 s l and solar zenith angle £50°) the dayglow will be the limiting 
factor in the detection of weak auroras. 


Introduction 

In a previous study, Germany et al. [1990] investigated the 
dependence of modeled far ultraviolet (FUV) atomic oxygen and 
N 2 Lyman-Birge-Hopfield (LBH) auroral emissions on likely 
uncertainties and anticipated changes in the neutral 
atmosphere. In this paper, we extend that study by examining 
the dependence of modeled FUV auroral emissions on the 
energy spectrum of the incident electron flux and by comparing 
the modeled FUV auroral emissions to anticipated underlying 
dayglow emissions. 

Germany ei al. [1990] modeled aurora] emissions (O I (1356 
A), LBH (1464 A), and LBH (1838 A)) to evaluate their utility 
in extracting fundamental information about the incident 
auroral electrons. As an example of this utility, a ratio of the 
column intensities formed from either of the shorter- 
wavelength emissions and the LBH (1838 A) emission (O I 
1356/LBH 1838 or LBH 1 464/LB H 1838) may be used to infer 
the characteristic energy of the incident auroral electrons. 
They also investigated the sensitivity of these emissions and 
their intensity ratios to changes in the modeled neutral 
atmosphere. The O I (1356 A) to LBH (1838 A) intensity ratio 
exhibits a greater dependence on the characteristic energy of 
the incident electrons than does the LBH (1838 A) to LBH 
(1464 A) intensity ratio. This is offset, however, by the 


denier for Space Plasma and Aeronomic Research, University of 
Alabama, Huntsville. 

*NASA Marshall Space Flight Center, Huntsville, Alabama. 
3 Department of Physics, Center for Space Plasma and Aeronomic 
Research, University of Alabama, Huntsville. 

4 Computer Science Department, Center for Space Plasma and 
Aeronomic Research, University of Alabama, Huntsville. 

Copyright 1994 by the American Geophysical Union. 

Paper number 93JA02357. 

0 1 48-0227/94 /93JA-02357I05. 00 


greater stability of the LBH emission ratio to changes in 
seasonal and solar cyclic activity. 

We have evaluated the sensitivity of selected FUV auroral 
emissions to a variety of energy distributions. Since different 
physical processes lead to different energy distributions, there 
is not a single parameterization which applies to all situations. 
Strickland et al. [1983] examined the dependence of FUV 
auroral emissions on the choice of either a Gaussian or a 
Maxwellian incident energy spectrum. Here we extend that 
study by modeling auroral emissions resulting from incident 
energy distributions of different widths and shapes to determine 
whether changes in energy spectra might obscure, or render 
ambiguous, our determination of the incident characteristic 
energy from FUV emission ratios. 

In addition, we investigate the relative intensities of 
modeled dayglow and auroral emissions as a function of solar 
activity levels and solar zenith angle. Since imaging in the far 
ultraviolet avoids problems with exceedingly bright visible 
dayside scattering, the detection limit is set by the brightness 
of the superimposed dayglow emission. 

Emission Models 

The auroral code used to compute the auroral intensities is a 
two -stream model described in an earlier paper [Rickards and 
Torr , 1990]. A two-stream model can serve well in cases such 
as the present which emphasize energy deposition and column 
intensities and in which detailed knowledge of pitch angle 
distributions is not crucial [. Stamnes , 1981]. As noted by 
Richards and Torr [1990], our auroral code gives excellent 
agreement with the more detailed multistream model of 
Strickland et al. [1983] and Daniell and Strickland [1986] 

The dayglow emissions are computed using our field line 
interhcmispheric plasma (FLIP) model optimized for airglow 
studies [Torr et al ., 1990]. Emission cross sections for N 2 LBH 
are those of Ajello and Shemansky [1985]. The O ( 5 S) 1356 
excitation cross section is that of Stone and Zipf[ 1974] scaled 
by the same factor of 3.1 by which their 1304-A emission 
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cross sections were reduced by Zip/ and Erdman [1985], 
Dissociative excitation of O z by electrons is a significant 
source of O I (7774 A) emission for hard auroral spectra 
[Strickland et at., 1989]. Since the same mechanism could also 
serve as an additional source of 1356-A emission, we modeled 
the effects of this additional source, using the cross sections of 
Wells et at. [1971] revised to reflect the lifetime determination 
°f Wells and Zipf [1974], We find that the increase in the 
auroral 1356-A column intensity with the additional excitation 
source is about 30% at 10 keV. The production of 1356-A 
dayglow by photoelectron dissociative excitation of 0 2 is 
considerably smaller, since photoelectrons are typically 
produced with energies below that of the peak e+0 2 cross 
section. For completeness we include the dissociative 
excitation source in the auroral 1356-A emissions modeled 
below. 

For the present study, both the auroral and the FLIP codes 
used the MSIS-86 neutral atmosphere [Hedin, 1987]. The 
model parameters differ slightly between the two studies 
reported below and will be given separately in each section. 
When necessary to model emissions for a range of activity 
levels, we vary the F \qj and Ap indices used in our emission 
models. The range of modeled activity levels is given in Table 
1 . 

We model three emissions as seen from a nadir-looking 
orbiting platform: O I (1356 A), LBH (1,1) (1464 A), and LBH 
(2,8) (1838 A). These emissions are chosen as representative 
of the principal FUV emissions available to space-based 
observations and as potential remote diagnostics. The promi- 
nent O I (1304 A) emission is not included because it is 
optically thick and undergoes multiple scattering. The 1356-A 
emission does undergo multiple scattering, but its efficiency is 
relatively small [Strickland and Anderson, 1983] and is ignored 
for the nadir-looking geometry considered here. In photo- 
metric imaging, the O I (1356 A) emission is not resolved from 
the LBH (3,0) (1354 A) emission which can be a significant 
contributor to a measured 1356-A signal (Figure 1). Hence we 
also model the LBH (1354 A) emission. Where appropriate, we 
present results for both the O I (1356 A) emission alone and for 
the summed O I (1356 A) and LBH (1354 A) emissions. In 
practice, we correct the signal at 1356 A for the LBH (3,0) band 
by measuring the longer- wavelength LBH bands separately. 

Dependence on Incident Energy Distribution 

We studied the dependence of FUV auroral emissions on 
selected energy distributions for the incident auroral electron 
flux by varying the energy spectrum and comparing resultant 
modeled emissions. In addition to varying the functional 
dependence of the energy distribution, we looked at the effects 
of varying its width and its low-energy behavior. 

We first modeled auroral emissions using either a Gaussian 
or Maxwellian energy distribution. Model parameters which 
might affect the modeled intensities, such as energy flux and 
average energy, were held constant throughout this study. 


TABLE 1. Input Paramete n 


Activity 


Minimum 

Moderate 

Maximum 

F 1&7 cm flux index 

75 

no 

200 

Ap magnetic activity index 

4 

20 

100 


Furthermore, care was taken to prevent high-energy truncation 
of the broader distributions due to either improper altitude or 
energy limits within the model. Emissions were modeled for a 
1 erg cm' 2 S'l incident energy flux; the moderate activity level 
from Table 1; geographic latitude = 60°; geographic longitude 
= 0°; solar apparent time = 0.0 hours, and day =173. 

After Daniell and Strickland [1986], we define the Gaussian 
energy distribution as 




Q 

it} 5 WE} 



( 1 ) 


and the Maxwellian distribution by 





( 2 ) 


where Q sets the incident energy flux. W sets the width of the 
Gaussian distribution, and E c is the characteristic energy of the 
incident distribution. The Gaussian scale parameter W is held 
constant at 0.25 for comparison with the Maxwellian 
distribution. The choice of this parameter is not important in 
itself, as we are not trying to match a particular auroral energy 
distribution at this point. The average energy <E> of an energy 
distribution f{E) is taken to be 

/ Ef(E) dE 

<E> =r- ( 3 ) 

I /(E) dE 

Figure 2 shows the two distributions for an average energy 
of 5 keV and a 1 erg cm' 2 s* 1 incident flux. The Maxwellian 
distribution is asymmetric, and its average energy will 
therefore not correspond to its characteristic energy (energy at 
which the distribution is maximum). Rather, the average 
energy of the Maxwellian distribution is equal to twice its 
characteristic energy. On the other hand, the Gaussian 
distribution is symmetric, and its average energy is thus equal 
to its characteristic energy. The average, not the 
characteristic, energy was held constant when comparing 
distributions. 

In Figure 3 we show the ratio of the auroral column 
intensities as calculated with the Gaussian and Maxwellian 
distributions. The longer-wavelength LBH (1838 A) emission 
is virtually unaffected (less than 5%) by the choice of energy 



Average Energy (keV) 

Fig. 1. OI 1356 and LBH 1354 A auroral emiirionj for a 1 erg on' 2 
r* incident energy flux and a moderate activity level as given in Table 1 . 
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Energy (keV) 

Fig. 2. Gaussian and Maxwellian energy distnbutions with a 5-keV 
average energy and a 1 erg on* 2 !* 1 incident energy flux. 




Average Energy (keV) 

Fig. 4. Same as Figure 2 except without 0 2 absorption. 



0 2 4 6 • 10 


Average Energy (keV) Energy (keV) 

Fig. 3. Ratio of integrated column brightne.iei calculated from a ener » r Attribution, with varying width.. Each 

Gau.iian energy distribution to those from a Maxwellian distribution. Attribution has an average energy of 5 keV and an energy flux of 1 erg 


distribution. Because of 0 2 absorption, the shorter- 
wavelength emissions exhibit significant variation with 
energy distribution. When the 0 2 absorption mechanism is 
disabled within the auroral code, we find that the LBH (1464 A) 
(and LBH (1354 A)) variations are entirely due to 0 2 absorption 
(Figure 4). The O I (1356 A) emission, however, is only 
partially dependent on the effects of O 2 absorption. 

The differences between modeled emissions from the 
Gaussian and Maxwellian distributions are due, in part, to the 
differences in the functional dependence of the two energy 
distributions on the incident energy, i.e., their shapes. In 
addition, however, the Maxwellian distribution is broader (in 
energy) than the Gaussian distribution. A wider energy 
distribution will contribute relatively greater numbers of both 
high- and low-energy electrons than a more narrow 
distribution. Since much of the energy dependence of FUV 
auroral emissions is due to absorption by 0 2 , some of the 
differences seen in Figure 3 can be expected to be due solely to 
the relative widths of the energy distributions. 

Accordingly, we selected a single distribution, a Gaussian, 
and modeled auroral emissions which differed only in the width 
of the distribution used. The width of the Gaussian distribution 
was varied by adjusting the width parameter W in equation (1) 
from 0.1 to 0.75. Sample distributions for an average energy 
of 5 keV are shown in Figure 5. Again, each distribution is 
constrained to have an energy flux of 1 erg cm* 2 i* 1 . Figure 6 




Average Energy (keV) 

Fig. 6. Ratio of integrated column brightnesses calculated from a 
Gaussian energy distribution with W=0.1 to those with W=0.75. 


shows column intensity ratios for W=0.75 and W=0.1. Note 
that the behavior is similar to that of Figure 3. 

We also studied auroral emissions using two arbitrary energy 
distributions. The intent of this comparison was to study the 
modeled emissions from energy distributions with quite 
different functional dependences and widths but with Fixed 
average energies. The first of these arbitrary distributions was 
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a step function which extended from 0 eV to a maximum energy 
equal to twice its average energy. The second distribution was a 
very narrow (VP=0.06) Gaussian distribution. Both distribu- 
tions, labeled “wide” and “narrow”, are shown in Figure 7 for an 
average energy of 5 keV and a 1 erg cm' 2 s* 1 incident flux. The 
choice of W for the narrow distribution results in a distribution 
whose width (full width at half maximum) is only 10% of its 
characteristic energy. With this combination of energy 
distributions, the wide distribution is about 20 times broader 
than the narrow distribution. Figure 8 shows the ratio of 
column emissions modeled with the two distributions. Despite 
the dissimilarities in energy distribution, the results shown are 
similar to those in Figures 3 and 6. Each of these studies 
implies that for a fixed average energy and energy flux, the 
width of the energy distribution is more significant than its 
shape in determining modeled FUV auroral column intensities. 

Previous work has pointed out the significance of the low- 
energy part of the energy spectrum in modeling visible 
emissions. For example, Meier et al. [1989] noted that 
measured energy distributions often have significantly more 
low-energy electrons than are provided by either the Gaussian 
or Maxwellian distributions given here by equations (1) and 
(2). This low-energy tail adds little to the total energy flux but 
should enhance F region emissions. The addition of the low- 
energy tail is not a significant emission source for the FUV 
emissions modeled here and results in changes of less than 
about 5%. 


Comparison to Auroral Dayglow 


The second task of this study was to compare the intensity of 
modeled auroral emissions with that of anticipated dayglow 
emissions as they would be seen from space for the features of 
interest here. Though FUV auroral images have served as 
remote diagnostics in a number of studies [e.g., Robinson et 
al., 1989. 1992], comparatively little work has been done for 
images seen against a full midday airglow background. We 
have recently developed a FUV imager for the Global Geospace 
Science mission (M. R. Ton et al.. An ultraviolet imager for 
the International So lar-Tenes trial Physics Mission, submitted 
to Reviews in Space Science , 1993) which incorporates 
ultraviolet filters both sufficiently narrow to allow separation 
of the FUV features of interest and with approximately 10 
orders of magnitude blocking throughout the visible and near 
infrared. We here assess the in-band competing elements of 
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Fig. 7. ‘‘Narrow'’ and “wide” energy distributions with a 5-keV average 
energy and a 1 erg cnr 2 r 1 incident energy flux. 


aurora and airglow by modeling FUV auroral emissions against 
a superimposed dayglow as a function of solar zenith angle and 
incident energy flux. 

Since the incident auroral flux is highly dynamic, individual 
measurements of incident energy fluxes may span a large range. 
For our study we desire an estimate of average energy fluxes 
coupled with an estimate of their anticipated frequency of 
occurrence. Such a study has been reported by Hardy et al. 
[1985] in which they analyzed data from two Defense 
Meteorological Satellite Program satellites as well as the 
Satellite Test Program P78-1 satellite. They reported average 
energy fluxes (keV cm- 2 s* sr 1 ) as a function of magnetic local 
time, magnetic latitude, and Kp for a period from September 
1977 to August 1980. We can use these results to estimate the 
average incident energy flux (ergs cm' 2 s 1 ) in specific parts of 
the auroral oval. The average energy flux reported by Hardy et 
al. [1985] for Kp £ 3 was less than about 10 ergs cnr 2 s*l in the 
midnight sector and less than about 3 ergs cm' 2 s* 1 in the noon 
sector. For the roughly 20% of analyzed data in which Kp 
exceeded 3, the average energy flux rose to about 25 ergs cnr 2 
s 1 in the midnight sector but stayed relatively constant in the 
noon sector. 

In a similar study, Spiro et al. [1982] analyzed data from the 
Low Energy Electron experiments on the Atmosphere Explorer 
C and D satellites from January 1974 to April 1976. They 
sorted their data by both Kp and the auroral electrojet index 
(AE) but presented detailed tabular results only for the AE 
index. They found that over 90% of their data occurred for 
values of AE less than 600. Within this selection, the average 
energy flux was less than about 6 ergs cm* 2 r 1 in the midnight 
sector and less than about 4 ergs cm -2 r 1 in the noon sector. 
The higher-activity cases showed an increase to about 9 ergs 
cm* 2 s* 1 for the midnight sector with no increase in the noon 
sector. This is essentially in agreement with Hardy et al. 
[1985], though the latter reported higher fluxes in the midnight 
sector. 

We have modeled incident auroral fluxes as Gaussian 
distributions (Ws0.25) with an incident energy flux of 1, 5, 
and 10 ergs cm* 2 S' 1 . Characteristic energies between 0.5 and 
10 keV were used for the incident auroral electrons. Local 
ionospheric conditions were modeled with the FLIP model 
using a flux tube characterized by an L shell of 5 for midsummer 



Average Energy (keV) 

Fig. 8. Ratio of integrated column brightnesses calculated from the 
narrow energy distribution shown in Pig. 7 to those from the wide 
distribution. 
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conditions (day=173). The solar zenith angle varied from 
about 50° at local noon to 83* near local midnight. To 
monitor the variability of the combined auroral/dayglo w 
intensity as a function of solar cyclic activity, the Fjo 7 and Ap 
input parameters were varied as shown in Table 1. The initial 
case modeled here (solar zenith angle of 51° at L=5, local 
summer noon) represents something of a worst case as auroras 
will typically lie to higher latitudes and hence larger solar 
zenith angles and less dayglow. The addition of the LBH (1354 
A) emission to the O I (1356 A) has little effect on the 
conclusions drawn from this part of the study. Therefore we 
omit the addition of the LBH (1354 A) in the remainder of the 
figures. 

Figure 9 shows a comparison of noon airglow and auroral 
emissions, where the latter have been computed for 1, 5, and 10 





Fig. 9. Integrated column brightnesses for auroral (solid) and dayglow 
(dashed) emissions with a solar zenith angle of 51°. All emissions are 
given for the range of activity levels defined in Table 1. Within each set, 
the largest values are obtained for maximum activity conditions. The 
auroral emissions are modeled for 1, 5, and 10 ergs cm* 2 ** 1 incident 
energy fluxes: (») O I 1356, (b) LBH U64, and (c) LBH 1838. 


ergs cm 2 s - i cases. The dayglow is strongly dependent on the 
selected activity level; the auroral is much less dependent. (For 
each wavelength, both the dayglow and the three auroral cases 
are modeled for three distinct activity levels.) 

For solar zenith angles of $51° the dayglow intensities 
exceed those of a 1 erg cm’ 2 r* aurora. Weak auroras with high 
characteristic energies would be the most difficult to detect at 
noon. Fortunately, dayside auroras tend to be caused by soft 
cusp particles. Figure 10 shows the ratio of 1356-A aurora to 
airglow for the 5 ergs cm" 2 s- 1 case for this same location as a 
function of solar zenith angle. As would be expected, the ratio 
is a strong function of solar zenith angle because of the airglow 
dependence on the latter. 

Discussion 

FUV auroral emissions have been modeled as a function of 
several different incident electron energy distributions. Our 
choice of distributions was not necessarily realistic; the 
distributions were chosen to illustrate extremes of shape and 
width. It is significant, in light of these choices, that a change 
of energy distribution affects the column emissions by less 
than 309b. Regardless of the energy distributions modeled, we 
see that the differences in resulting emissions are governed 
primarily by the average energy of the distributions and the 
relative widths of the distributions. Therefore, much of the 
behavior seen in Figures 3, 6, and 8 can be discussed solely in 
terms of the relative widths of the energy distributions. This is 
principally due to the predominate influence of O2 absorption 
in determining the energy dependence of the FUV auroral 
emission. 

The relative width of the energy distribution affects the 
modeled emissions via either a low-energy or a high-energy 
mechanism, with the distinction between the two being at an 
average incident energy of about 10 keV. Below 10 keV, FUV 
auroral emissions generated by the broader distribution 
experience a greater loss at the low-altitude edge of the 
emission altitude profile and produce subsequently reduced 
column intensities. For average energies above 10 keV, 0 2 
absorption dominates, and only those emissions due to the 
lower-energy part of the distribution are seen. (The emissions 
from the higher-energy electrons are totally absorbed.) 
Longer-wavelength emissions, such as LBH (1838 A), do not 
experience 0 2 absorption and are relatively unaffected by the 
choice of incident energy distribution. The O I (1356 A) 



Characteristic Energy (keV) 

Fig. 10. Ratio of modeled dayglow tad auroral 1356-A emu nous ai a 
function of solar zenith angle for a 5 ergs cm* 2 s* 1 energy flux and a 
moderate activity level as given m Table 1. 
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emission exhibits an additional dependence due its altitude 
distribution. A broad energy distribution will result in reduced 
1356- A production, relative to a narrow distribution, because 
of increased collisional competition with other neutral targets. 

As noted above, the magnitude of the changes in column 
intensities with choice of energy distribution is less than about 
30%. The dependence of modeled FUV auroral emissions on 
energy distribution is therefore much less than their 
dependence on characteristic (or average) energy. The char- 
acteristic energy extracted from FUV color ratios should 
therefore be relatively unambiguous, since such determinations 
will not be dependent on an exact knowledge of the incident 
energy spectrum. 

Typical aurora have also been modeled with concurrent 
dayglow emissions as a function of solar zenith angle, solar 
activity levels, total electron energy flux, and incident 
characteristic energy. The ratio of auroral to dayglow 
emissions is a strong function of solar zenith angle because of 
the dayglow dependence on solar zenith angle. For the FUV 
emissions of interest here, the dayglow is the strongest 
competing background and under certain conditions (energy 
flux <1 erg cm" 2 s' 1 and solar zenith angle ^50°) will be the 
limiting factor in the detection of weak auroras. 
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We have studied the behavior of modeled auroral 
Hall and Pedersen conductivities as a function of 
auroral electron flux and atmospheric 
composition. We show that the column 
integrated Hall and Pedersen conductivities may 
be determined from a single measurement of 
two FUV auroral LBH emissions: LBH 1464 A and 
LBH 1838 A. The dependence of the determined 
conductivities on incident characteristic energy, 
total energy flux, changes in solar activity 
levels, and arbitrary perturbations of the 
individual neutral components (O, O 2 , N 2 ) have 
been sutdied. The conductivity determinations 
are stable with respect to changes in levels of 
solar activity and against changes in the neutral 
atmospheric composition. The incident energy 
flux is also provided by the LBH 1838 A 
emission. Thus the two auroral LBH emissions 
are sufficient to determine the column 
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Abstract 


The purpose of this paper is to examine the viability of using LBH emission ratios to infer 
auroral conductances and to quantify the strengths and weaknesses of this technique. We 
show that column-integrated Hall and Pedersen conductances may be determined from a 
single remote measurement of a pair of auroral LBH emissions, one in the region of strong 
O2 absorption (1464 A) and one lying outside of this region (1838 A). The dependence of 
the determined conductivities on incident average energy, total energy flux, and changes in 
solar and magnetic activity levels is examined. We show that, for energies above 5 keV, 
auroral conductances may be scaled by the square root of the incident energy flux with 
errors less than about 20%. For lower energies, however, the scaling may deviate 
significantly from a power of 0.5. We provide appropriate scaling factors as a function of 
average energy. We also note that the choice of either a Gaussian or Maxwellian 
distribution can significantly affect low-energy (<1 keV) conductance calculations. Finally, 
we quantify the conditions under which the mean energy of a Gaussian energy distribution 
may differ significantly from the characteristic energy. 



introduction 


Knowledge of ionospheric conductivities is central to an understanding of the coupled 
nature of the ionosphere, thermosphere, and magnetosphere. This is reflected in models of 
ionospheric electrodynamics and convection requiring conductance distributions [e.g., 
Richmond and Kamide, 1988; Rasmussen and Schunk, 1987; Senior and Blanc, 1984; 
Spiro and Wolf, 1984; Roble et al., 1982]. Considerable effort has been directed toward 
modeling ionospheric conductances from EUV solar illumination, nocturnal sources, and 
auroral precipitation [see, for example, Moen and Brekke, 1993; Brekke and Hall, 1988; 
Rasmussen et al., 1988; and the reviews by Brekke and Moen, 1993 and Reiff, 1984]. 
Conductivities have been inferred from ground-based radar observations by a number of 
investigators, including Robinson et al. [1987,1985]; Robinson and Vondrak [1984]; and 
Vickrey et al. [1981, 1982]. Satellite data is commonly used to determine conductances 
[e.g., Watermann et al., 1993; Fuller-Rowell and Evans, 1987; Vondrak and Robinson, 
1985; Wallis and Budzinski, 1981], with many studies featuring satellite overflights of 
ground-based facilities. 

A particular difficulty in modeling conductances from auroral precipitation is proper 
accounting of the high temporal and spatial variability in the auroral zone, while at the same 
time mapping its global morphology. This difficulty can potentially be overcome by using 
images of the auroral region taken from space. For example, auroral images from 
Dynamics Explorer have been used to infer ionospheric conductances [Robinson et al., 
1989; Kamide et al., 1986]. In addition, Robinson et al. [1992] and Basu et al. [1993] 
used FUV emissions measured by the AIRS instrument on the Polar BEAR satellite [Meng 
and Huffman, 1987; Schenkel and Ogorzalek, 1987] to infer E region electron density 
profiles. 

A key to remote determination of auroral conductances by optical methods is the ability 
to characterize the incident auroral flux. Since the resultant column intensities are relatively 



insensitive to choice of energy distribution [e.g., Germany etal., 1994] the principal task is 
to estimate the energy flux and characteristic energy. Rees and Luckey [1974] showed that 
measurements of the 6300 A emission of atomic oxygen and the 4278 A emission of 
ionized molecular nitrogen could be used to infer the characteristic energy of the incident 
auroral particle flux. In this technique, the energy dependence is provided by the 
collisional deactivation of the 6300 A emission at lower altitudes. Strickland et al. [1983] 
demonstrated the same capabilities using FUV auroral emissions. In particular, they 
proposed the use of the O I (1356 A) and N 2 Lyman-Birge-Hopficld (LBH) emissions as 
remote diagnostics of the auroral process. Instead of quenching mechanisms, the FUV 
technique relied on absorption of O 2 at lower altitudes to supply an energy dependence. 
The potential of using FUV emissions from the Viking satellite as remote measures of 
characteristic energy was examined by Steele et al. [1992]. They found only a weak 
dependence on energy, but they were only able to examine energies below 1.4 keV. 
Robinson et al. [1992] compared electron density measurements from the Sondrestrom 
radar with electron densities inferred from O I (1356 A) and N 2 + (3914 A) emissions 
observed by AIRS and found good agreement between the measured and inferred profiles. 
Since their method depends on a visible emission, it cannot be used for sunlit regions 
because of jhe intense Rayleigh scattering. To allow sunlit observations, Strickland et al. 
[1983] proposed the use of two FUV emissions, O I (1356 A) and LBH (1273 A), as 
diagnostics. 

Germany et al. [1990] examined the stability of FUV emission ratios formed from 
atomic oxygen and LBH emissions as a function of likely changes in the neutral 
atmosphere. They found that emission ratios formed from the O I (1356 A) emission 
varied significantly with seasonal and solar cyclic variations. On the other hand, emission 
ratios employing long- and short-wavelength LBH emissions were much less sensitive to 
changes in the neutral atmosphere. Furthermore, blending of the O I (1356 A) and LBH 
(1354 A) emissions represents a serious observational difficulty. Consequently, the LBH 
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emission ratios are potentially more valuable as remote diagnostics of the incident auroral 
flux than similar ratios using the O I ( 1 356 A) emission. Germany et al. [ 1994] examined 
the use of FUV emission ratios as remote diagnostics as a function of incident energy 
distribution. They found that the intensity of modeled FUV emissions are governed 
principally by the mean energy and relative width of the energy distribution. Column- 
integrated intensities are relatively insensitive (<30% ) to changes in the actual functional 
form, i.e., the shape, of the energy distribution. Consequently, the LBH emission ratio 
can be used as remote diagnostics of auroral flux without a priori knowledge of the energy 
distribution. Robinson et al. [1987] noted similar results for auroral conductances, i.e, 

they are primarily dependent on the mean energy and energy flux, and less sensitive to the 
energy distribution. 

This paper is a continuation of our previous examinations of the utility of LBH 
emission ratios as remote diagnostics of auroral processes. Here, we examine the 
relationship between FUV LBH auroral emissions and ionospheric conductances due to 
auroral precipitation. To our knowledge, this is the first tune the LBH emission ratio has 
been used to infer auroral conductances. It is the purpose of this paper to examine the 
viability of using LBH emission ratios to infer auroral conductances and to quantify the 
strengths and weaknesses of this technique. We demonstrate that determinations of 
average energy and energy flux from a single measurement of two LBH emissions may be 
used to estimate column-integrated ionospheric conductances. In addition, we examine 
common assumptions regarding the energy flux scaling and the average value of a 
Gaussian energy distribution, and examine the conditions of their validity. 

Models 

Conductances are directly proportional to the corresponding ionospheric electron 
densities. While FUV emissions may be used to infer characteristics of the auroral 
processes, they must be coupled with additional information to provide the electron 
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densities needed for the conductance determination. This information may be provided by 
ground-based radar measurements, by in situ measurements, or by numerical models. We 
calculate conductances using the densities, temperatures, and magnetic field strengths 
provided by the field line interhemispheric plasma (FLIP) model [Torr et al., 1990]. The 
FLIP model solves the coupled time-dependent equations for energy, momentum, 
continuity, and photoelectron and auroral transport along a magnetic field line extending 
between the northern and southern hemispheres. The complete chemistry in the model is 
presented by Torr et al. [1990] and Torr [1985]. 

Incident auroral fluxes are modeled with a two-stream electron energy-loss model 
[Richards and Torr, 1990] and input to the FLIP model. The two-stream auroral model 
gives excellent agreement with the more detailed multistream calculations of Strickland et 
al. [1983] and Daniell and Strickland [1986] The auroral code is optimized through the 
use of variable energy and altitude grids to calculate the degraded energy spectrum of the 
incident electron flux and the subsequent emission of a number of prominent auroral 
emissions, including O I (1356 A) and the Lyman-Birge-Hopfield band system of N2. The 
incident auroral particle flux is assumed to be electrons only; precipitating ion fluxes are not 
modeled. Both the auroral and FLIP models use the MSIS-86 thermospheric model of 
Hedin [1987]. 

For this study, the incident auroral electron fluxes are modeled by Maxwellian energy 
distributions with average energies between 500 eV and 20 keV and energy fluxes between 
1 and 10 erg cm* 2 s*L Additionally, for comparison, calculations using Gaussian energy 
distributions are made for the same range of average energy and energy flux. Local 
ionospheric conditions are modeled for a flux tube characterized by an L- shell of 5 and 
magnetic longitude of 180° for midwinter conditions (day=356) and near local midnight 
(local time=1.4 hours; solar zenith anglc=127°). As a consequence, the ionospheric 
conductances are assumed to be solely of auroral origin, i.e., there is no EUV contribution. 
In geographic coordinates, the flux tube is located at (74N, 11 IE) in the northern 



hemisphere and (51S, 1 1 IE) in the southern hemisphere; the results presented are for the 
northern location. We have modeled different activity cases by adjusting the F\ 0.7 and Ap 


model inputs. Maximum conditions are represented by Fio. 7=200 and Ap=lOO; moderate 
conditions by Fio. 7=110 and Ap=20\ and minimum conditions by F 10 . 7=75 and Ap= 4. 
Results are presented for the moderate case except where noted. 

The Pedersen volume conductivity, op, is a function of the charged particle mobility in 
the direction of the local electric field and perpendicular to the local magnetic field. The 
Hall conductivity, 0 //, is dependent on the mobilities in the direction perpendicular to both 

the E- and fi-fields. 


<*P = f | »ekp + X n i k P J 

(la) 

oh = J- | n e kfi - X n i k l H 

(lb) 

The mobility coefficients are functions of the ratios of the collision frequency v and the 

gyrofrequency C2 . 
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Our model includes three collisional interactions: ion-neutral (subscript ik in equations 
2), electron-neutral ( ek ), and electron-ion (ei). The electron-neutral collision frequencies 



arc from Schunk and Nagy [1978]. The electron-ion and ion-neutral collision frequencies 
are from Schunk and Nagy [1980] except for the 0-0 collision frequency, which is taken 
from Salah [1993]. In practice, the conductivities below 250 km are dominated by the ion- 
neutral collision frequencies. The electron collision frequencies make no significant 
contributions to the conductivities, except above 250 km where the electron-ion collision 
frequency contributes slightly. For completeness, we include all three collision frequencies 
in the results discussed below. 

Ionospheric conductivities are dependent on ionization due to both solar EUV 
illumination and incident particle flux. The combined FLIP and auroral models provide 
ionization from both sources. Brekke and Hall [1988] noted potential difficulties from not 
separating auroral and non-auroral ionization sources in conductance calculations. As 
noted above, we have chosen winter midnight conditions to minimize the solar EUV 
ionization sources as we are primarily interested in this study in the correlations of 
conductance with emissions from auroral particle fluxes. However, there are still nighttime 
sources from scattered light and starlight. To ensure we are dealing only with auroral 
ionization sources, we have subtracted underlying, i.e., non-auroral , ionization sources 
from the total ion production. Similar corrections are performed for the auroral emissions. 
This correction for non-auroral sources is significant principally for incident electron fluxes 
with characteristic energies below 1 keV. 

Moen and Brekke [1990], using incoherent scatter radar measurements, derived quiet- 
time conductances due to solar EUV ionization and concluded that ion composition 
variations have little effect during solar minimum conditions. They noted that the height- 
dependent conductivities can be well represented without the explicit ion summations used 
in Eqs. (1). The present study, however, focuses only on auroral ionization sources. 
Low-energy auroral conductivities may be quite sensitive to ion composition. Therefore, 
the ion densities used are those of the major ion species (NO+, C>2 + , N2 + , O-, and N+). 
Local charge neutrality is assumed to hold; hence the electron density is taken to be the sum 
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of the ion densities. Only the three prominent neutrals are considered: O 2 , N 2 , and O. The 
ions are assumed to have a common temperature. 

The column integrated conductances are the sum of the volume conductivities over the 
column height 

Zpm = j Gpji dz ( 3 ) 

The integration in Eq. (3) spans the entire range of available data (90 to 500 km). Column 
integrated conductances are a function of the average energy of the incident particle flux 
(Figure 1). The figure also shows conductances from Vickrey et al. [1981] and 
Rasmussen et al. [1988] for comparison. The three studies essentially agree for the 
Pedersen conductance with some differences for the Hall conductance. Rasmussen et al. 
[1988] noted that the choice of thermospheric model can significantly influence the modeled 
conductances. Vickrey et al. [1981] used the 1000 K thermosphere of Banks and Kockarts 
[1973] while Rasmussen et al. [1988] used the MSIS-83 model [Hedin, 1983]. We use 
the MSIS-86 model [Hedin, 1987]. Brekke and Hall [1988] used a simplified analysis to 
show that the Hall-to-Pedersen conductance ratio is directly proportional to k — v# / n*. 
Therefore, differences in the Hall profiles may represent differences in choice and 
calculation of collision frequencies, in addition to differences in thermospheric models. In 
the discussions below, we use the term conductivity for the local volume conductivity and 
the term conductance will be reserved for the column integrated conductivities as given in 
Eq. (3). The units of volume conductivity are mho/m (siemens/m in SI units); those of 
conductance are mho (siemens). 

Dependence on energy flux 

The volume conductivities are dependent on the electron and ion densities which, in 
turn, are functions of the total incident energy flux. The continuity equation for the ion 
density, N,-, is given by 





(4) 
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where Qi is the ion production rate, L, is the corresponding loss rate, and the final terra 
represents transport processes. As used here, the subscript i identifies total production and 
loss rates for all ion states. The principal losses are dissociative recombination, 

NO + + e =N +0 

N 2 + e =N + N (5) 

Oj + c —O + 0 

and ion-atom interchange, 

0 * + 0 2 = 02+0 

2 ( 6 ) 

o+ + n 2 =no + + n 

[Torr, 1985]. For illustration, we assume the recombination reactions of Eq. (5) proceed 
with a common recombination rate, a. Similarly, we assume the interchange reactions of 
Eq. (6) can be described by a common reaction rate, k. Under these assumptions, the loss 
rate, L„ is given by 

Li = o*e]([Atf +fNO + ]+[0j]) + k[0 + ]{[N2\+[0 2 }) (7) 

where square brackets are used to denote number density of a given species. 

If we ignore the transport term in Eq. (4) and assume equilibrium conditions, the 
production and loss rates can be equated. In the FI region at about 150 km, the 
dissociative recombination rate is the rate limiting process in Eq. (7). Furthermore, the 
sum of ion species in Eq. (7) is equal to the electron density at this altitude, so the 
production rate can be written as 

<2 = o(e] 2 


(8) 



where we have dropped the subscript i in recognition of the fact that the total ion and 
electron densities are equal. In the F2 region between 200 and 400 km, the ion-atom 
interchange reactions of Eq. (6) are the rate limiting process. Furthermore, the O is the 

dominant ion species and can be taken as equal to the electron density. Under these 
conditions the production rate is 


Q = k[e]([N 2 M0 2 ]}. (9) 

In both cases, the electron density can be solved for in terms of the production rate, 

[e] = cQ r ( 10 ) 

where c is a proportionality constant and /ranges from 0.5 at lower altitudes to 1.0 at 
higher altitudes. Since the ionization rate Q is proportional to the incident energy flux and 
the conductances are linear functions of the ion densities, the col umn conductances will 
scale with the energy flux, i.e.. 


Epjt - cEjf ( 11 ) 

where Ef is the incident energy flux and c is, again, a proportionality constant Since the 
conductivities typically peak in the lower thermosphere, it is usually assumed that /=0.5. 
For example, Moen and Brekke [1993] and Robinson and Vondrak [1984] noted that 
column-integrated conductances due to only EUV ionization vary with the square root of 
the 10.7-cm solar flux, which is assumed to be proportional to the incident energy flux. A 
value for /of 0.5 is strictly valid, however, only below about 150 km where dissociative 
recombination is the important loss rate. 

Since we are discussing only auroral ionization sources, the ionization profiles peak at 
different altitudes depending on the average energy of the incident particle flux. Hence a 
single scaling will not necessarily be appropriate for all average energies. This is especially 
significant for relatively narrow energy distributions, such as the Gaussian used for 
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comparison below. Broader distributions, e.g. Maxwellian, will produce ionization across 
a greater altitude range, resulting in smaller variations in /with average energy. In Figure 
2 the column-integrated Pedersen conductivity is shown as a function of the incident energy 
flux for several average energies. Results for both Gaussian (dashed line) and Maxwellian 
(solid line) distributions are shown. In Table 1 we give the exponent /from Eq. (1 1) as a 
function of the average energy of an incident electron flux. The exponent was obtained 
from a fit to column conductivities modeled with incident energy fluxes ranging from 1 to 5 
erg cm-2 s-1. For energies below about 2.5 keV /exceeds 0.5 by more than 20% for the 
Pedersen conductances. 

Dependence on Energy distribution 

The choice of modeled incident energy distribution can also be significant for low-energy 
fluxes, as shown in Figure 3. In the figure we show Hall conductivity profiles modeled 
for Gaussian and Maxwellian distributions with an average energy of 500 eV and an 
incident energy flux of 1 erg cm* 2 s* 1 . We define the Gaussian energy distribution as 



and the Maxwellian distribution as 

/max(£) = exp -(-£-)] (13) 

where Q sets the incident energy flux, W sets the width of the Gaussian distribution, and 
E c is the characteristic energy of the incident distribution. These are the same distributions 
used by Strickland et al. [ 1983] except that we have replaced the Gaussian scale parameter 
W of Strickland et al. [1983] with WE C in Equation (12). As used here, W is a 
dimensionless parameter. For comparison with the Maxwellian distribution, we set the 
Gaussian scale parameter W=0.25. The choice of this parameter is not important in itself 
as we arc not trying to match a particular auroral energy distribution at this point. This 



value was chosen to match that of Strickland et al. [1983]. The average energy <£> of an 
energy distribution f(E) is taken to be 
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( £ ) = - 2 • (14) 

mdE 

Robinson et al. [1987] discussed the dependence of calculated conductivities on the 
choice of energy distribution and addressed potential confusion arising from the use of 
different energy distributions. They pointed out that the average energy of the asymmetric 
Maxwellian distribution is twice its characteristic energy. The Gaussian energy distribution 
used in Figure 3, however, is symmetric and the average and characteristic energies are the 
same. In an appendix, we examine conditions under which the average and characteristic 
energies of a Gaussian energy distribution may differ significantly. For a given average 
energy, the high-energy tail of the Maxwellian results in more ionization at lower altitudes 
than from the Gaussian used here with IV=0.25 (Figure 4). The additional ionization 
between 110 km and 180 km results in increased Hall conductivity over the same altitude 
range. The difference in conductivity for the two energy distributions is not significant for 
incident fluxes with average energies above about 5 keV. 

The choice of energy distribution also affects modeled FUV auroral emissions, 
although to a much smaller degree. Germany et al. [1994] examined FUV auroral 
emissions as a function of energy distribution and found relatively small (<30% ) 
dependence for column-integrated emissions. The observed differences were governed 
principally by line-of-sight loss mechanisms, i.e., absorption by O2 in the Schumann- 
Runge absorption continuum (1300 - 1750 A). Consequently, the differences are largest 
for average energies of 10 keV. 




Robinson et al. [1987] obtained good agreement between conductances calculated with 
energy spectra measured by the Hilat satellite and those based on a Maxwellian energy 
distribution. We use the Maxwellian distribution given in Eq. (13) without attempting to 
include more detailed modifications, e.g., a low- or high-energy tail [Meier et al., 1989; 
Strickland et al., 1993]. To facilitate comparison with other energy distributions, we 
present results as a function of average energy (twice the Maxwellian characteristic energy). 
If the maximum modeled energy in our auroral code is set too low there will be significant 
truncation of the high-energy tail of the Maxwellian distribution. We avoid this potential 
problem by setting the maximum modeled energy to be 100 keV, well above the largest 
average energy (20 keV) modeled here. 

The curve labeled “nighttime” in Figures 3 and 4 is due to non-auroral ionization 
sources included in the FLIP model, principally from resonantly-scattered sunlight, 
airglow, and starlight [Strobel et al., 1980]. The nighttime ion production rates are 
typically much smaller than those from auroral ion production. For incident fluxes with 
low average energies, however, the electrons do not penetrate into the E-region, leaving the 
nighttime photoionization as the dominant ion source in the lower thermosphere. The 
contribution to the conductivities from nighttime photoionization can be significant for low- 
energy auroras. As noted above, the nighttime ionization sources are subtracted from the 
results shown below. 

CONDuenvmES from Auroral emission Ratios 

The principal requirement for remote determination of both the Pedersen and Hall 
conductivities is a remote diagnostic of the characteristic energy and of the energy flux of 
the incident auroral flux. Such diagnostics can be provided from ratios of selected LBH 
bands [Germany et al., 1990]. In Figure 5 we show the column-integrated intensity of two 
LBH auroral emissions (LBH (1464 A) and LBH (1838 A)) as a function of characteristic 
energy for a 1 erg cm- 2 s* 1 incident energy flux. The LBH emission ratio (Figure 6) is a 
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clear function of characteristic energy and, unlike similar ratios using 01 emissions, is 
virtually independent of seasonal and solar cyclic changes in neutral composition. The ratio 
thus represents a particularly useful tool for the determination of the characteristic energy of 
precipitating electrons. 

The utility of these emission ratios is based on the absorption of FUV auroral emissions 
via the O 2 Schumann-Runge absorption continuum (1300 - 1750 A). Because the LBH 
emissions are excited solely by electron impact the LBH column-intensity for a nightside 
aurora is a measure of the total energy flux, provided there are no significant loss 
mechanisms, e.g., O 2 absorption, in the viewing column. Longer-wavelength LBH 
emissions, such as LBH 1838, can thus be used as remote diagnostics of energy flux while 
shorter-wavelength, energy-dependent emissions (LBH 1464) may be used in conjunction 
with the longer-wavelength emissions to monitor average energy. The ratio of LBH 1838 
and LBH 1464 serves as an indicator of characteristic energy that is independent of the 
incident energy flux. 

With the use of the LBH emissions as a remote diagnostic of characteristic energy and 
energy flux, the integrated conductivities may be determined. In essence, we replace the 
average energy scale in Figure 1 with the LBH emission ratio and then scale by the incident 
energy flux as determined from the long-wavelength LBH measurement This is shown in 
Figure 7 in which the conductances are plotted as a function of the LBH emission ratio. 
(The individual points are modeled for average energies of 0.5, 1, 2.5, 3.5, 5, 7.5, 10, 
12.5, 15, 17.5, and 20 keV.) In Figure 8 the dependence of conductivities derived from 
the LBH emission ratio is shown as a function of activity. The activity levels were varied 
by changing the F 10.7 and Ap indices, as discussed above. Since the LBH emission ratio 
is relatively insensitive to these changes, most of the variation seen in Figure 8 is due to 
changes in the ionization and temperature profiles. 

It is at the low energies that the present method of inferring the conductivities is least 
sensitive because the LBH ratio is relatively insensitive to characteristic energies below 2 
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kcV. In addition, the conductivities are more sensitive to solar activity changes changes in 
energy distribution at low energies. Furthermore, though they were removed for this 
study, nighttime ionization sources can exceed ion production in soft auroras. Clearly care 
must be taken when modeling lower energy auroras. These difficulties may conceivably 
be addressed by employing other remote diagnostics with greater sensitivity to low 
energies, such as the O I (1356 A) to LBH (1838 A) ratio. Unfortunately, using the O I 
(1356 A) line to obtain increased sensitivity to low energies has the disadvantage of 
increased sensitivity to thermospheric perturbations. In addition, working with the 1356- A 
emission is complicated by additional difficulties, e.g., blending with the LBH (1354 A) 
emission. 

SUMMARY 

In this study the calculation of ionospheric conductivities is linked to remote 
observables (FUV LBH emissions) which could be observed by a space-based, nadir- 
looking instrument with selected bandpasses in the far ultraviolet. The study is limited to 
auroral particle sources and examination of EUV ionization sources is deferred to a future 
study. A single observation of a pair of LBH emissions, one with wavelength less than 
about 1750 A and the other at longer wavelengths, is sufficient to characterize both the 
incident energy flux and its characteristic energy. This characterization is relatively 
insensitive to changes in the ionosphere and thermosphere due to solar cyclic variations. 
Just two parameters, energy flux and characteristic energy, convey sufficient information to 
characterize the incident particle flux and, in turn, to infer the resultant ionospheric 
conductivities. The LBH emission ratio is least sensitive to changes in energy for average 
energies below 1 keV. In addition, the low-energy conductance calculations are 
particularly dependent on the choice of energy spectrum. We also examine the square-root 
scaling rule with incident energy flux usually applied to conductances. For our modeled 
cases, auroras with average energies greater than about 5 keV can be scaled by the square 
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root of the energy flux with errors less than about 20%. However, below 5 keV the 
scaling may deviate significantly from a power of 0.5, dependent on the modeled incident 
energy distribution. 
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APPENDIX 

It is common to state that the mean energy of a Gaussian energy distribution is its 
central, or characteristic, energy. Since energy is constrained to be non-negative, however, 
the distribution may be truncated at zero energy. If so, the distribution will be assymetric 
and the mean energy will be larger than the characteristic energy. For example, in Fig. Al, 
the curve labelled with w=0.75 is truncated at zero. Therefore its mean energy is 
somewhat larger than the 5 keV mean shared by the other distributions. The discussion 
below details the conditions under which the mean energy of a Gaussian distribution may 
be taken as equal to its characteristic energy. 

The Gaussian distribution is defined in terms of two parameters, the central value jJ. and 
a width factor <r. 

/*(£;/i,o)=Ae-(vf (Al) 

The distribution in Eq. (Al) is the same as that given in Eq. (12) with /i=£ c and <r=wE c . 
The width factor is related to the full width at half maximum (T) by setting fg(£+/72; fx, o) 
= 0.5 f g (u; ii, a). 

r = lilnla (A2) 

The mean energy is the weighted average across the distribution. 



where the limits of integration reflect the constraint that energy be non-negative. 

The integral in the denominate* - is evaluated first, using a transformation to the variable 


z 
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I 2 =J f g (E)dE 
Jo 

= A<rJ e* 2 dz 


(A4) 


where 


G 

dE = crdz 


and y is the ratio of the central energy and the width parameter. 




The integral is evaluated in terms of the error function, 

I 2 = A<x< 


-J^ e* 2 dz + j e * 2 dz 


= Aa& -{erf(>^ + 1 } 


(A5) 


(A6) 


(A7) 


where 


crf(z)s-£r[ e'\* 
and erf(-z) = - erf(z) 


(A8) 


The integral in the numerator of Eq. (A3) is now evaluated, again using the transformation 
from Eq. (A5). 
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■1 


I,= EfJE) dE 


l 


= Ac r| (<rz + n)e- z dz 
-T 


-Ad 2 I ze zl dz 

J-y 




(A9) 


Note that the integral is given in terms of I2, the integral in Eq. (A7). To complete the 
evaluation, a second transformation is used. 


y = z 2 
dy = 2 zdz 


(A10) 


Ii 



e-ydy + nh 


ASL e -r> + 

2r 


I 2 


(All) 


In the last step of Eq. (All), the definition of yis used to introduce a common factor of fj. 
in both terms. 

The mean energy can now be written, using Eqs. (A3), (A7), and (A1 1), as 

{£) = |i- = /iF(» (A 12) 

where F()) is given by 


F(tf = 1 + 


sit . 

Vny{erf(>)+ 1} 


(A 13) 


Figure A2 shows F(>) for values of /between 0.1 and 10. As can be seen, for values 
of /greater than 1.5 the low energy truncation of the distribution is negligible and the mean 



energy can taken to be the characteristic energy. For values of y below 1.5, the truncation 
is significant and the approximation that <E> = Ec is not valid. Since y is the ratio of the 
central energy to the width of the distribution, we see that the approximation of the mean 
energy by the characteristic energy fails for energy distributions that are very wide 
compared with the central energy. 

In this work, we set a = wfi. The requirement that ybe greater than 1.5 therefore 
translates to the condition that w be less than 0.667. For our work, w = 0.25 and the mean 
energy may thus be safely approximated by the characteristic energy. 
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Fig. 1. Column-integrated Hall and Pedersen conductances for an incident energy flux of 1 
erg cm- 2 s _1 and a Maxwellian energy distribution. Dashed line, Vickrey et al„ 
[1981]; solid line, Rasmussen etal., [1988]; circles, this work. 

Fig. 2. Column-integrated Pedersen conductances as a function of incident energy flux and 
average energy. Dashed line, Gaussian; solid line, Maxwellian. 

Fig. 3. Hall conductivity due to an incident auroral flux with an average energy of 500 cV 
and an incident energy flux of 1 erg cm-2 s* 1 . Conductivities are modeled for both 
a Gaussian and a Maxwellian energy distribution. The curve labelled “Nighttime” 
is due to non-auroral sources. 

Fig. 4. Electron density profiles corresponding to the Hall conductivity profile in Figure 3. 

Fig. 5. Column-integrated intensity of LBH (1838 A, 1464 A) for a 1 erg cm* 2 s-l incident 
energy flux and a Maxwellian energy distribution. 

Fig. 6. Ratio of auroral LBH (1838 A) to LBH (1464 A) for a 1 erg cm* 2 s-1 incident 
energy flux. Three cases are shown: Fio.7=75, Ap=4; F io.7=l 10, Ap=20; 
Fio.7=200, j4p=100. 

Fig. 7. Column-integrated Hall and Pedersen conductances as a function of auroral LBH 
emission ratio. The individual points are modeled for average energies of 0.5, 1, 
2.5, 3.5, 5, 7.5, 10, 12.5, 15, 17.5, and 20 keV and a 1 erg cm- 2 s* 1 incident 
energy flux. 

Fig. 8. Column-integrated Hall and Pedersen conductances as a function of auroral LBH 
emission ratio for a 1 erg cm- 2 s- 1 incident energy flux. Three cases are shown: 
Fiq.7=75, Ap=4 (dot); Fio.7=100, Ap= 20 (line); Fio.7=200, Ap = 100 (dash). 
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Figure A 1 . Gaussian energy distributions with fixed characteristic energy and varying 
widths. The plots are labelled with the dimensionless parameter w = a/fi = l/y. 

Figure A2. F(# as a function of y. F(1.5) = 1.0202 
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TABLE 1. Power-law dependence of conductances with energy flux. 


y (Maxwellian) /(Gaussian) 


Average Energy (keV) 

Pedersen 

Hall 

Pedersen 

Hall 

0.5 

0.74 

0.60 

0.80 

0.62 

1.0 

0.66 

0.55 

0.76 

0.65 

2.5 

0.58 

0.54 

0.60 

0.57 

5.0 

0.56 

0.53 

0.57 

0.55 

10.0 

0.54 

0.52 

0.54 

0.53 

15.0 

0.54 

0.52 

0.54 

0.52 
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Fig. 1. Column- integrated Hall and Pedersen conductances for an 
incident energy flux of 1 erg cm* 2 s* 1 and a Maxwellian 
energy distribution. Dashed line, Vickrey etal., [1981]; 
solid line, Rasmussen et al., [1988]; circles, this work. 
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Fig. 2. Column-integrated Pedersen conductances as a function of 
incident energy flux and average energy. Dashed line, 
Gaussian; solid line, Maxwellian. 
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Fig. 3. Hall conductivity due to an incident auroral flux with an 
average energy of 500 eV and an incident energy flux of 1 
erg cm-2 s* 1 . Conductivities are modeled for both a 
Gaussian and a Maxwellian energy distribution. The curve 
labelled “Nighttime” is due to non-auroral sources. 
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Fig. 4. Electron density profiles corresponding to the Hall 
conductivity profile in Figure 3. 
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Fig. 5. Column-integrated intensity of LBH (1838 A, 1464 A) for a 
1 erg cm-2 s-1 incident energy flux and a Maxwellian energy 
distribution. 
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Fig. 6. Ratio of auroral LBH (1838 A) to LBH (1464 A) for a 1 erg 
cm* 2 s' 1 incident energy flux. Three cases are shown: 
^ 10 . 7 = 75 , Ap=4; F\qj= 100, Ap=20) F io.7=200, Ap= 100. 
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Fig. 7. Column-integrated Hall and Pedersen conductances as a 
function of auroral LBH emission ratio. The individual 
points are modeled for average energies of 0.5, 1, 2.5, 3.5, 
5, 7.5, 10, 12.5, 15, 17.5, and 20 keV and a 1 erg cm-2 s*l 
incident energy flux. 
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Fig. 8. Column-integrated Hall and Pedersen conductances as a 

function of auroral LBH emission ratio for a 1 erg cm* 2 s* 1 
incident energy flux. Three cases are shown: Fio.7=75, 
Ap=A (dot); F] o. 7 = 100, Ap=20 (line); Fio.7=2CO, -4/7=100 
(dash). 
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Figure A 1 . Gaussian energy distributions with fixed characteristic 
energy and varying widths. The plots are labelled with the 

dimensionless parameter w = o//i = 1/y. 




Figure A2. F(# as a function of y. F(1.5) = 1.0202 



